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ABSTRACT 

We report the UV luminosity function (LF) of Lyman break galaxies at z ~ 5 derived 
from a deep and wide survey using the prime focus camera of the 8.2m Subaru telescope 
(Suprime-Cam). Target fields consist of two blank regions of the sky, namely, the 
region including the Hubble Deep Field-North and the J0053+1234 region, and the 
' total effective surveyed area is 1290 arcmin 2 . Applications of carefully determined 

colour selection criteria in V — I c and I c — z' yield a detection of 853 z ~ 5 candidates 
^ ■ with z' AB < 26.5 mag. The UVLF at z ~ 5 based on this sample shows no significant 

Qh' change in the number density of bright (L > L* =3 ) LBGs from that at z ~ 3, while 

there is a significant decline in the LF's faint end with increasing lookback time. 
£^ ■ This result means that the evolution of the number densities is differential with UV 

^ \ luminosity: the number density of UV luminous objects remains almost constant from 

C3 • z ~ 5 to 3 (the cosmic age is about 1.2 Gyr to 2.1 Gyr) while the number density 

of fainter objects gradually increases with cosmic time. This trend becomes apparent 
thanks to the small uncertainties in number densities both in the bright and faint 
parts of LFs at different epochs that are made possible by the deep and wide surveys 
we use. We discuss the origins of this differential evolution of the UVLF along the 
cosmic time and suggest that our observational findings are consistent with the biased 
galaxy evolution scenario: a galaxy population hosted by massive dark haloes starts 
active star formation preferentially at early cosmic time, while less massive galaxies 
increase their number density later. We also calculated the UV luminosity density 
by integrating the UVLF and at z ~ 5 found it to be 38. 8j^'^% of that at z ~ 3 
for the luminosity range L > 0.1L* =3 . By combining our results with those from the 
literature, we find that the cosmic UV luminosity density marks its peak at z = 2-3 
and then slowly declines toward higher redshift. 

Key words: cosmology: observations - galaxies: high-redshift - galaxies: luminosity 
function - galaxies: evolution - galaxies: formation. 



1 INTRODUCTION 

The galaxy luminosity function (LF) is one of the funda- 
mental quantities frequently used to investigate the prop- 

* Based on data collected at Subaru Telescope and partly ob- 
tained from the SMOKA science archive at Astronomical Data 
Analysis Center, which are operated by the National Astronomi- 
cal Observatory of Japan, 
f E-mail: iwata@oao.nao.ac.jp (II) 



erties of galaxy populations and explore the evolutionary 
process of galaxies. It represents the number densities of 
objects with specific luminosity at specific redshift ranges. 
Recent development of large imaging and spectroscopic sur- 
veys improved the accuracy of LFs in the local universe in 
optical and near-inf r ared wavelengths dNorberg et al.l 120021 : 
iBlanton et ail 120031 ; iKochanek etafl l200ll ). Goodaualitv 
LFs at UV wavelengths have also been obtained by com- 
binations of balloon or spaceborne missions with ground- 
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based survey data (|Sullivan et al.ll2000l ; IWvder et al.ll2005h . 
For LFs at earlier cosmic times, in principle a complete 
and volume-limited sample of galaxies at the target red- 
shift is required. Great efforts have been made by spectro- 
scopic surveys to obtain large g alaxy samples capable of de- 
riving LFs up to z < 2 (e.g., iLillv et all 1 19951 : iLin et all 
Il999l ; lllbert et al.ll2005l ). However, a construction of a large 
spectroscopic sample of galaxies at higher redshift, z > 4, 
would require enormous observing time due to the faint- 
ness of the target objects and, consequently, it remains a 
big issue to be addressed by future observing facilities. In- 
stead of determining redshifts of large numbers of objects, 
we are able to derive a luminosity function statistically by 
using a sample of galaxies with redshift estimates based 
on their spectral energy distributions traced by multi-band 
photometry. In this way, extensive multi-wavelength sur- 
veys conducted by the Hubble Space Telescope (HST) and 
large ground-based telescopes have enabled us to measure 
LFs at different epochs up to z < 2-3 (e.g.. [ Sawicki et aF 



19971 ; iDrorv et al 



20031; iKashikawa et al.ll2003T 
200:!: I'ozzotti all 12003! : IWolf et all 120031 

2002 



2( 



Chen et al 



Dahlen ct al 
Poli et alJ i 2003 



2005) and earlier (e.g.jLanzetta et al..,~, ,, ... , . ... . , _ , 

Gabasch et"aH l2004al ; iGabasch et alj 120061 ; iPaltani et al 



2006), and to discuss the evolution of galaxies through the 



comparison of LFs at different cosmic times. 

One of the most powerful observational methods to 
probe galaxies at high redshift (z > 2) is the detection 
of the spectral discontinuity due to the redshifted Lyman 
limit (912A) or Lyman a (caused by absorption by inter- 
galactic hydrogen gas) through multi-wavelength broad- 
band imaging. This method is called the Lyman break 
method, and si nce the pioneering work b y Steidel and 
his co -workers (jSteidel and Hamilton! 1 19921 ; ISteidel et al.l 
Il996al ). large number of Lyman brea k galaxies (LBGs) 
at z ~ 2-4 have been detected (e.g . . iMadau et al. 
Steidel et al.ll2003l; ISteidel et al J 12004 ; iFoucaud et al 
Sawicki and Thompsonll2005l ). 



1996 



2003 



There are several obstacles in searching for higher red- 
shift (2 > 4) LBGs. First, the larger distance makes the 
apparent brightness of target objects dimmer so that, for 
the same luminosity, an object at z = 5 is about 1 magni- 
tude fainter than that at 2 = 3. Next, the sensitivity of CCD 
chips declines at A > 8000 A. And finally, for ground-based 
instruments, the increase of background sky brightness due 
to the night sky lines emitted from the Earth's atmosphere 
makes the S/N much worse than in the bluer wavelength 
ranges. Despite these difficulties, it is important to push to 
these earlier epochs in the history of galaxy formation. 

In Iwata et al. (2003; hereafter Il03h we created a sam- 
ple of « 300 LBGs with J C (AB) < 26.0. It was the first 
large sample of LBGs at z ~ 5, and was made possi- 
ble by the unique combination of the large mirror aper- 
ture of the 8.2m Subaru telescope and the wide field 
of view of its prime focus camera named Suprime-Cam 
ijMivazaki et al.l l2002h . It was also in this work that the 
UVLF of LBGs at z ~ 5 was derived statistically for 
the first time. We found no significant change in the 
bright luminosity range (Muv J; —21 mag) between z — 



5 and z = 3. After II03I several studies have success- 



fully created LBG samp l es at 2 ~ 5 (|Le 



inert and Bremerl 



20031 ; lOuchi et a' . l l2004al ; [Giavalisco et al.l I l2004bl ; iLee et al.l 
2006). Attempts to obtain statistically meaningful sam- 



ples of LBGs at 2 > 6 have also been fruitful 

fe.g.. Bouwens et alj 20031; Bouwens and Illingeorthl 120061; 



iBouwens et all 120061; Istanwav et all 120031 ; I Dickinson et al.l 
|2004| ; iBunker et alj|2006h . As we accumulate knowledge of 
the physical propertie s of high-redshift (2 > 2) galaxies 
such t heir as sizes (e.g..[steidel et"al"1ll996al ; lGiavalisco et all 



19961; Ravindranath et alT 20061). ste l lar p o pulations (e.g 



Sawicki and Yedll998l; IShaplev etail l200ll; iPapovich etTal] 



erties (e.g., Giavalisco et al. Il998l; 


Adelberger et al. 


» 1 — 1- 
1998; 


Ouchi et al.l 2004b; Hamana et alj 


2004; Kashikawa ct al. 


2006), we improve our chances of constructing a compre- 



hensive view of galaxy evolutionary processes in the high 
redshift universe and eventually through cosmic time. 

One serious limitation of past ground-based LBG sur- 
veys is that only luminous objects have been studied. The 
statistical nature of galaxies fainter than t he typical lumi- 
nosity (L* in the LF parameterized with the ISchechterlll976l 
form) has been left unexplored, although these numerous 
faint galaxies are estimated to dominate the rest-frame UV 
luminosity density. Thus, a comparison between the evolu- 
tion of luminous objects and that of faint objects has not 

been made. 

The Keck Deep Fields (KDF) (|Sawicki and Thompson] 

120051) made a very deep imaging survey reaching 7?.(AB) ~ 
27 mag and constructed samples of LBGs at z ~2 to 4 in- 
cluding objects far fainter than L* . They find that the UVLF 
at z ~ 4 shows a significantly smaller number density in 
the fainter part than that at z ~ 3, while in the bright 
end no difference between z ~ 4 and 3 has been found. 
Such differential evolution of the UVLF was not detected in 
the preceding LBG work with brighter limiting magnitudes 
(|Steidel et alj|l999h . The KDF indicated that deep surveys 
reaching fainter than L* (with sufficient survey area) are 
required to closely inspect the nature of high-z galaxy evo- 
lution. It is an interesting subject to examine whether such 
differential evolution started at earlier epoch, z > 5. For 
that subject the construction of deep sample of LBGs down 
to 1 magnitude fainter than L* at z > 5 is indispensable. 

Here we report the UVLF of LBGs at 2 ~ 5 based on our 
updated sample that builds on the work of 103. The survey 
now consists of two independent fields and has an effective 
survey area of 1290 arcmin 2 after masking bright objects. 
The limiting magnitudes of the sample LBGs are z = 26.5 
mag for one field (the field including the Hubble Deep Field 
- North) and 25.5 mag for the other (J0053+1234). The 
sample limiting magnitude is now ~ 0.5 mag deeper than 
in our previous work (103), which was 26.0 mag in I c . The 
survey area is more than two times larger than our previous 
survey, which enables us to examine the degree of field-to- 
field variance and to estimate how representative are our 
measured properties of 2 ~ 5 galaxies. Through the com- 
parison of the z ~ 5 UVLF obtained with this deeper and 
wider sample with the deep U VLFs at z ~ 3 and 4 by 
ISawicki and Thompson! (|2006al ). we attempt to disclose the 
evolution in the UVLFs of LBGs. 

In section 2 we describe the observations for the two 
target fields. Procedures of data processing, catalog con- 
struction and sample selection are described in section 3. 
In section 4 we present our 2 ~ 5 UVLFs and also compare 
them with previous results. In section 5 we make a com- 
parison of our 2 ~ 5 UVLF with those at 2 = 4 and 3 in 
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ISawicki and Thompsonl (|2006ah and then discuss the differ- 
ential evolution of UVLF. Then in section 6 we derive the 
UV luminosity density at 2 ~ 5 and discuss the evolution 
of the UV luminosity density through the cosmic time. Our 
findings are summarized in section 7. Throughout this pa- 
per we adopt a cosmology with the parameters fijvf = 0.3, 
fi A = 0.7, H = 70 km/s/Mpc. Magnitudes are on the AB 
system. 



2 OBSERVATIONS 

We chose two target fields for our survey of LBGs at 
2 ~ 5. One is i n the area of the H ubble Deep Field 

19961), 
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North (HD F-N; IWilliams et al 



J0053+1234 (|Cohen et al.lll996l ; ICohen et ailll999al ). Both 



fields are target fields of existing extens i ve redshift surveys 
(e.g.. ICohen et all l2000t iDawson et ail l200ll ; IWirth et all 
12004 ; fcowie et alj|2004 ICohen et al.lll999bh . and there are 
many galaxies with spectroscopic identifications at interme- 
diate redshifts (2 < 4). Such information is valuable to de- 
fine reliable colour selection criteria for LBGs at our target 
redshift range, z = 5. 

We use the V, I c and 2' filters to select our 2 ~ 5 
LBGs. In Fig. [T]we show the transmissions of these filters 
as a function of wavelength. In the figure we also indicate 
model SEDs of star-forming galaxies with and without dust 
attenuation that we generated using version 2 of the popula- 
tion synthesis code PEGASE l|Fioc and Rocca-Volmerangel 
Il997li A consta nt star formation rate with the Salpeter IMF 
l Salpetedfl955ft is assumed, and the age from the onset of 
star formation is 100 Myr. Dust attenuation we show here 
is E(B — V) = 0.4 using the attenuation law proposed by 
ICalzetti et ail (120001 ) and absorption by intergal actic neutral 
hydro gen is applied using the analytic formula in llnoue et ail 
(2005). The redshift we assumed here is 4.8, which is an av- 
erage redshift of our sample LBGs estimated with simula- 
tions (see section 4.1). In this figure the model SEDs are 
normalized at 10,000 A in the observer frame. The VI c z' 
filter combination efficiently detects the break at Lycv for 
star- forming galaxies at 4.3 < 2 < 5.3, as verified by simula- 
tions (described in section 4.1) and f ollow-up spectroscopy 
l|Ando et ail |2004| ; Undo et ail 120071 ). In addition to our 
VI c z' filters, Fig. QJalso shows t he t'-b and filter (dotted line) 
and we note that Ouchi et ail (|2004al ) use the combination 
of V, i' and 2' filters to select their LBGs at 2 ~ 5. The 
i'-band filter is shifted slightly shortward compared to the 
7 c -band filter used in our observations. 
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Figure 1. Effective transmissions (including quantum efficiency 
of the CCD and atmospheric absorption) for the three filters used 
in our observations (solid lines with blue (V), green (J e ) and 
red(2 ; )). The transmission of i' filter used in lOuchi et alj ll2004al ) 
is also shown as a dotted line. Model spectra of star-forming galax- 
ies at z = 4.8 are shown as black lines. The solid line represents 
the case without dust attenuation and the dashed line is for the 
case with E(B - V) = 0.4. These SEDs are normalized at 10,000 
A. See text for details of the models. Upper-side ticks indicate 
rest-frame wavelengths for an object at z = 4.8. 



from each other in several aspects. First, the central field po- 
sitions are ~3' apart, and in our observations the position 
angle was set to 10° (east of north), while in the observa- 
tions by the UH group position angles were 0°or 180°. The 
position and the position angle in our observation were ad- 
justed so that the HDF-N taken by the HST/WFPC2 was 
located at the centre of one of Suprime-Cam's CCD chips. 
The dithering lengths are also different. In our observations 
dithering scales are only 10"-15", which were chosen so that 
the area with deepest integration time became as large as 
possible. We could choose such small dither lengths because 
we concentrated on the detection of galaxies at 2 ~ 5, which 
sh ould be apparently small (< 2"). In the UH observations 
bv lCapak et ail (|2004l ) the dither scales were around l'and 
the effective exposure time within the observing field is fairly 
close to uniform. In this paper the area common to these two 
different sets of observations is used, since that area achieves 
the deep image depths in all of three observing bands (V, 
I c and 2'). In Table [1] we summarize observations for the 
HDF-N region. 



2.1 The Hubble Deep Field - North 

We made deep imaging observations of the field that sur- 
rounds the HDF-N using the Suprime-Cam attached to the 
8.2m Sub aru Telescope in February 2001. The results pre- 
sented in Il03l are based on these data alone. For the new 
results we show in the present paper we also used the imag- 
ing data taken by the University of Hawaii (UH) group 
l|Capak et alj|2004h using the same instrument and the same 
filters during February 28 - April 23, 2001. We retrieved 
these UH data from the SMOKA data archive jBaba et ail 
12002ft . 

Properties of these two sets of observations are different 



2.2 J0053+1234 

In order to increase the survey area and examine the effect of 
field-to-field variance, we carried out the imaging of another 
blank field. Specifically, we imaged the J0053+1234 region, 
which is one of th e target fields of the Caltech Faint Galax y 
Redshift Survey l|Cohen et ailll999al ; ICohen et ailll999bh . 
Using extensive Keck spectroscopy, these authors identified 
163 galaxies in the redshift range 0.173 < 2 < 1.44 within 
an 2' x 7/3 field. Steidel et al. (1999, 2003) also selected 
this region as one of the target fields of their searches for 
LBGs at 2 ~ 3 and 4, and they identified 19 galaxies with 
spectroscopic redshifts at 3 < 2 < 5. 
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Our imaging observations of this field using Subaru / 
Suprime-Cam have been executed during the summer sea- 
sons from 2002 to 2004. The central coordinate of the point- 
ing is 00 h 53 m 23.0 3 , +12°33'56"(J2000.0). In Table [2] we 
summarize information of our observations and final lim- 
iting magnitudes in V, I c and z'. Although the point spread 
function of the final images are smaller than that in our 
images of the HDF-N region, the total limiting magnitudes 
achieved for this field are 0.3-0.5 mag shallower than those 
for the HDF-N region since the effective net exposure times 
are shorter due to the poor weather conditions we experi- 
enced. Nevertheless, these data more than double the total 
survey area and so are very valuable for populating the rel- 
atively rare, bright part of the sample and for examining 
field-to-field variations. 



3 DATA REDUCTION AND SELECTION OF 
LBG CANDIDATES 

3.1 Image Processing 

The basic image processing was done using nekosoft - 
the software de veloped for the data reduction of mosaiced 
CCD cameras (|Yaei et all 120021 ') — as well as with IRAF. 
Bias was subtracted using overscan regions. Next, flat-field 
frames were generated by combining median-normalized ob- 
ject frames after masking objects. Then, after division by 
the flat-field frame and cosmic-ray removal (using CRUTIL 
in IRAF), a distortion correction was made for every frame 
with correction factors which were determined for each CCD 
position. 

Image mosaicing in each band was made by identify- 
ing 40-180 stars common in object frames. The correction 
for differences in sensitivity and fluctuations of point source 
flux was also made based on these star data. In the HDF-N 
region, the data taken by us have small dither scales and so 
we could not adopt a mosaicing procedure usually executed 
in nekosoft, which determines relative positions of frames by 
common stars in different CCD chips. Thus, we first derived 
relative positions of frames taken by the UH group and gen- 
erated a mosaiced image using the UH data alone. We then 
calculated the positions of the frames taken by us relative 
to that mosaiced image. The final mosaiced image was then 
made using both the UH and our data. Among the 10 CCD 
chips of the Suprime-Cam, one chip was not working during 
part of the observations of the HDF-N region. We also ex- 
cluded data of another CCD chip in our data because there 
is not enough overlapping area between the positions of this 
chip and the UH data to allow us to determine reliable rela- 
tive positions. In the J0053+1234 region data from all CCD 
chips were used. The effective survey areas (after masking 
bright objects) are 508.5 arcmin 2 (the HDF-N region) and 
781.4 arcmin 2 (the J0053+1234 region). The FWHM values 
of point sources in the mosaiced images are about l.l"for the 
HDF-N region and 0.9"for the J0053+1234 region (Tables[Q 
and [2]). 

3.2 Astrometry 

The conversions from pixel coordinates in mosaiced CCD 
images to equatorial coordinates were made using the 



USNO-B1 catalog jMonet et all 120031 ). In total 1,000 to 
1,300 stars in the catalog were identified in the mosaiced 
image in each band, and polynomial coefficients up to third 
order were calculated. Internal position errors in coordinate 
conversions are less than 0.2"rms all over the images. Be- 
cause these accuracies are almost identical to that of the 
USNO-B1 catalog itself, we expect that positional uncer- 
tainty due to the errors in the corrections for optical image 
distortions and those in the mosaicing process to be smaller 
than 0.2". 

3.3 Photometric Catalog 

The determination of ph otom etric zero points for the HDF- 
N region is described in Il03l . Briefly, we det ermined photo - 
metric zero points based on imaging data of lLandohJ (|l992T j 
standard stars for V and I c bands, while for z'-band the zero 
point was determined through the distribution of I c — z 
colours of field stars and galaxies exposed on our images. 
Systematic errors in the determination of zero points are 
estimated to be 0.02, 0.08 and 0.1 mag for V, I c and z'- 
bands, respectively. The photometric zero points of the final 
mosaiced image s we re determined by comparison with the 
images used in 103. For the J0053+1234 region, the same 
procedure was followed for the V and I c bands, while for 
z'-band images of spectrophotometric standard stars taken 
during the same observing nights were used. Systematic er- 
rors are estimated to be 0.04 mag for U-band and 0.03 mag 
for I c and z'-bands. Corrections for Galactic ext i nction were 
made based on the dust map bv lSchlegel et al.l (|1998l ). The 
amount of reddening for the HDF-N region is assumed to 
be E(B - V) = 0.012, and for the J0053+1234 region it is 
E{B-V) = 0.066. 

Object det ection and photometry w ere made by us- 
ing SExtractor (Be rtin and Arnoutsll 19961 ) ver. 2.3 in single- 
image mode. A Gaussian filter with FWHM=3 pixels was 
used, and at least 4 contiguous pixels with counts higher 
than 1.3 a above background counts were required to be de- 
tected as an object. The cross-identifications of objects in 
our three bands were made in equatorial coordinates. Ob- 
jects were included in the object catalog only when they 
were within the area with sufficient image depth, determined 
by checking the exposure maps which take account of the 
number of exposures and variations in sensitivity between 
individual CCD chips. This procedure guarantees that the 
image quality is fairly close to uniform and that the vari- 
ance of photometric errors does not affect the uniformity of 
the galaxy sample. For U-band dropout LBGs we require 
objects to be detected both in I c and z'-bands. We used 
MAG_AUT0 as a total magnitude in all bands, and derived 
V — I c and Ic — z' colours using 1.6"diameter aperture mag- 
nitudes. Since object centroids are determined separately in 
each band, there is a possibility that the positions of aper- 
tures in different bands have offsets (e.g., star- forming re- 
gions are prominent in 7 c -band and their distribution is ir- 
regular). We feel that such problems are rare, since FWHM 
sizes of LBG candidates are mostly smaller than the aper- 
ture size used to derive colours (1.6"). We also find through 
visual inspection that morphologies of U-dropout objects in 
7 c -band and in z'-band are similar in most cases. The pho- 
tometric catalogs of both fields were constructed based on 
the z'-band. For objects undetected by SExtractor in the 
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l^-band, l/-band aperture photometry was made at the ob- 
jects' positions in the z'-band. In these cases, if the flux 
density within 1.6"diameter exceeded the 3-cr sky fluctu- 
ation then the measured aperture magnitude was used to 
calculate V ~ I c colour. Otherwise the 3-cr upper limit was 
adopted. 

3.4 Selection of LBG Candidates 

The colour criteria we adopted to se lect LBG candidates at 
z ~ 5 are the same as those used in Il03l . They were defined 
to efficiently select star-forming galaxies at z ~ 5 without 
heavy contamination by objects from lower redshifts, and 
were chosen by considering the colours of both model galax- 
ies and real galaxies within our images which have spectro- 
scopic redshifts identified by previous surveys available at 
that time (see below for the cross-identifications of galax- 
ies with spectroscopic redshifts with our object catalogues). 
Our colour criteria are expressed by the following equations 
in the AB magnitude system: 

V - la > 1.55 
and 

V-I c > 7.0(I C - z) + 0.15. 

In Fig. [2] we show our z ~ 5 colour selection criteria 
along with redshift colour tracks of several types of galax- 
ies. The colour tracks of model star-forming galaxies are 
calculated based on the same SEDs as used in Fig. [T] (see 
section 2). Colour tracks for spiral (Sbc and Scd) and el- 
liptical gal axies are calcula t ed us ing template SEDs of local 
galaxies by I Coleman et all (|l980h . 

In Fig. [3] we show colours of detected objects within a 
magnitude range 23.0 < z < 24.5. The numbers of objects 
which satisfy the colour criteria are 617 and 236 for the 
HDF-N region (z' < 26.5 mag) and the J0053+1234 region 
(z' < 25.5 mag), respectively. 

Among the HDF-N field LBG candidates in ll03l that are 
within the region common to the new, deeper and narrower 
imaging data, all of the objects with I c ^ 24.5 mag and 82% 
of objects with I c ^ 26.0 mag are included in our updated 
source catalog. Th e nu mber of objects that were listed as 
LBG candidates in ll03l but are outside the colour criteria in 
the new catalog is 20-30% larger than the estimated number 
of low-redshift contaminants (see section 4.1 and Table [5]). 
This fact is not surprising, because the signal-to-noise ratios 
of our previous data are worse than the HDF-N region data 
presented in this paper. 

Both of our target fields have been observed with 
other intensive spectroscopic surveys. We identified objects 
with spectroscopic redshifts in our images and plotted their 
colours in Fig. El For the HDF -N region we used the result s 
of two surveys l|Cohen et all |2000| and IWirth et ai1l2004h 
The n umbers of ob jects identified are 644 for I Cohen et al.l 
(l2000h and 1.803 forlWirth et all J2004h (most of the obj ects 
in lCohen et al.ll2000l are included in lWirth et~alll2004h . In 
Fig. EJa) we only show objects with I c > 22.0 mag and 

V > 22.0 mag, because photometry of objects brighter than 
these limits would be unreliable due to saturation or detec- 
tor non-linearity. Out of 1237 objects with spectroscopic red- 
shifts smaller than 4.5, only two fall within our z ~ 5 colour 
selection criteria. Both of them have z'-band magnitude 



I 
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Figure 2. A V — I c and I c — z' two colour diagram. The 
shaded area stands for the colour selection criteria adopted for 
the selection of Lyman break galaxies at z ~ 5. Lines are loci 
of colour changes for model galaxies along redshift. Solid blue: 
a star-forming galaxy without dust attenuation. Solid orange: a 
star-forming galaxy with dust attenuation of E(B — V) = 0.4. 
Solid circles are plotted for z 4.1 with an interval of 0.1. Sym- 
bols at z = 4.5 and z = 5. are enlarged. Colour tracks calculated 
using template SEDs by IColeman et al.l |l980j) are also shown. 
Cyan (dot-dashed): Sbc, light blue (dotted): Scd, red (dashed): 
elliptical. Symbols are plotted with a redshift interval of 0.5. Star 
symbols indic ate the colours of AO— M9 stars calculated based on 
the library bv lPickles! Jl99Sh . 



» 22 and their redshifts are 0.32 (ID=7066 bv lWirth et~al] 
|2004 ) and 0.64 (ID=3249). For the J0053+1234 region, we 
show objects with published spectroscopic redshift smaller 
than 4.5 with magenta squares in Fig.|3jb). One sensible ex- 
ception to this is CDFb-G5 at z = 4.486, which satisfies our 
colour selection criteria (V — I c = 2.70, I c — z' = 0.29) and 
is shown with a filled green circle as one of the spectroscopi- 
cally confirmed LBGs at z ~ 5 (see section 3.5). In addition 
to this object, two objects at z < 4.5 are within the colour 
selection cri teria for z ~ 5 LBGs . One of them is object id 
DOK 163 inlCohen et all l|l999al ). The redshift provided in 
ICohen et all l)l999ah is z = 0.580 and is tagged with 'absorp- 
tion only; faint; id uncertain'. It has a blue colour in I c — z' 
(0.07) in our photo metry. Another object is CDFa-GD12 in 
ISteidel et al.l (|l999h . The object is selected with their z ~ 4 
LBG criteria (using G, 1Z, I filters) and the spectroscopic 
redshift is z = 3.469. Although the /-band magnitude listed 
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m ISteidel et all (|l999h is similar to our z'-band magnitude 
(z ' = 24.76 mag ) , the distance between the coordinate listed 
m ISteidel et afT (|l999h and ours are fairly large (1.74"), un- 
like other objects (mostly < 1"). 

In Il03l we defined the colour criteria for selection of 
z ~ 5 LBGs using the colours of galaxies with spectroscopic 
and photometric redshifts in the HDF-N region, as well as 
model colour tracks. Since the publication of 103 the number 
of low-redshift galaxies with spectroscopic identifications in 
the HDF-N region has greatly increased, but, as described 
above, the cross-identification of these spectroscopic samples 
with our LBG catalog shows that the number of low-redshift 
galaxies that contaminate our colour selection criteria is very 
small. 

Since the magnitude range of our LBG candidates is 
fainter than most of these spectroscopically identified low- 
redshift objects, the objects identified as contaminants can- 
not be used to estimate the fraction of contaminants in our 
sample. Consequently, we need resampling tests to estimate 
the contamination fraction as a function of magnitude, as 
described in section 4.1. 

Our HDF-N region contains the GOODS-N area, where 
deep HST/ ACS and Spitzer/IRAC a nd MIPS imaging have 
been made (jGiavalisco et al] l2004al ). The ACS images are 
in the -B435, Veoe, I775 and zg 50 filters and using Veoe, hi5 
and Z85o, a U-dropout L BG sample has been constructed 
l|Ravindranath et al.ll2006l ). However, it is not straightfor- 
ward to compare our sample with those based on the ACS 
images. Since the filters used are different, especially Veoe 
which has a very broad bandpass (width is ~ 1500A), the 
colour tracks of the same galaxy SEDs with GOODS/ACS 
filters in the two-colour diagram are quite different from 
ours. 

In Fig. [4] the sky distributions of LBG candidates 
in our two fields are shown. There is a clear fluctua- 
tion in surface density of LBG candidates; some areas 
have many LBGs while others don't, and these irregu- 
larity in the surface density is not caused by the masks 
used to avoid effect of foreground bright objects such as 
galactic stars and nearby galaxies. Fairly strong cluster- 
ing signals h ave been found with past wide-field LBG 
searches (e.g., 



Giavalisco et al.ll 19981 ; [Adelberger etafll 19981; 
20031 : lOuchi et all l2004bl ; lHamana et all 



Fouc aud et al. 

20061 ). These figures also tell us that in small-area surveys 



such as the HDF-N, the Hubble Ultra Deep Field and even 
GOODS, the effect of field-to-field variance would be much 
more serious than for our data. Detailed analyses of clus- 
tering properties of our LBG samples will be presented in a 
forthcoming paper. 



3.5 Spectroscopic Identifications of Galaxies at 



We carried out optical spectroscopic observations us- 
ing a multi-obje ct spectroscopy mode of the FOCAS 

l|Kashikawa et all I2002T ) attached to the Subaru telescope. 
Until now, 8 objects in the HDF-N region and 4 ob- 
jects in the J0053+1234 region have been identified to 
be in the redshi f t range between 2 = 4.2 and 5.2 
i|Ando et all |2004| ; lAndo et all 120071 ; including one AGN 
GOODS J123647.96+620941.7). No object with smaller red- 
shift has been found in our observations so far, although 



more than half of the observed spectra have too low a 
signal-to-noise ratio to determine their redshifts. Four ob- 
jects which have slightly redder I c — z colours than the 
colour sele ction criteria have been identified as Galactic Mi- 
type stars (|Ando et al.ll2004l ). 

In addition, there are 8 objects in the HDF-N region 
and 2 objects in J0053+1234 which have been spectroscop- 
ically identified to be at 4.5 < z < 5.5 by other researchers. 
Among them, 6 objects are selected as U-dropout candidates 
in our images. Tables [3] and [3] provide summaries of the 
cross-identifications of these objects with previously pub- 
lished spectroscopic redshifts. In Fig. [3] the objects that 
have been selected as U-dropout candidates in our images 
and that have spectroscopic confirmations are shown with 
filled green circles (or upper arrows for objects not de- 
tected in U-band). There are four objects in the HDF-N 
region which have been previously reported to have spec- 
troscopic redshifts larger than 4.5 and are detected in our 
images but do not satisfy our colour selection criteria (see 
table |3]). Their positions in the V — I c and I c — z' di- 
agram are shown as open circles in Fig. [3] The first of 
these, HDF-N 4-439.0 at z = 4.54, has V - I c = 1.48 
and is slightly off from our colour selection window. An- 
other object, F36219— 1516, which has been reported to be 
at 2 = 4.89 bv lDawson et al.l (|200ll ). has a very blue V — I c 
colour in our image (V — I c — 0.67). Since the redshift deter- 
mination was made with solely a single emission line (identi- 
fied as Lya), we suggest that this may be a misidentification 
of a lower redshift object. F36376— 1453 is another object 
with a single emission line in th e spectrum wh i ch wa s iden- 
tified as Lya at z = 4.886 by iDawson et all (|200ll ). This 
object is very bright (z' — 21.97) and its colours in both 
V -Ic and I c - 2' are red (V - I c = 3.25, I c - z' = 0.78 and 
this object is outside of plotted area of Fig. [3]). It has a round 
shape in our images in I c and z'-bands (FWHM~ 1.0"). The 
fourth object, GOODS J12 3721.03+621502.1, was identified 
as an object at z = 4.76 bv lCowie et al] (|2004l ) (no. 1173 in 
their table 1). This object also has red colours in V — I c and 
I c — z . No information about t he quality of the spectrum is 
provided in lCowie et al] (|2004l ). On the other hand, one ob- 
ject i n the J0053+1234 r egion which is at z = 4.486 (CDFb- 
G5 in ISteidel et al] [l999) and is within our selection window 
(V — Ic = 2.70, I c — z = 0.29; see tabled, and it is shown as 
an filled green circle in Fig.[3jb). In summary, spectroscopic 
observations for our LBG candidates still cover only a small 
fraction of the sample, but the results obtained so far are 
encouraging. Both the existence of objects at z ~ 5 within 
the selection window and the small number of contaminants 
described in section 3.4 strongly support the validity of our 
selection criteria. 

ICapak et al.l l|2004l ) analysed multicolour data includ- 
ing some of the same Suprime-Cam data that we now re- 
processed and included in the present work to construct our 
multi-band object catalog in the HDF-N region. They also 
made a selection of U-dropout galaxies, deri ved their num- 
ber counts, a nd compared them w ith those in ll03l The num- 
ber counts in ICapak et ail (|2004T ) are just 13-21% of those 
in this paper. They claimed that our sample may be con- 
taminated by 2 ~ 1 galaxies since our selection criteria are 
redder than theirs in I c — z' colour. However , their argu- 
ments are incorrect. The colour criteria used in lCapak et al] 
(2004) are V -I c ^ 2.4 and V - I c ^ 7(I C -z')- 0.2 and the 
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Figure 3. V — I c and I c — z' two colour diagrams, (a): data points are for objects in the HDF-N region. Black circles indicate colours 
of objects with 23.0 < z' ^ 24.5. For objects not detected in V-band, arrows are used to show the 3-cr lower limit s of V — I c colours . 
Red and blue squares r epresent galaxies wi th spectroscopic redshifts (z < 4.0) published in 



Ked and blue squares represent galaxies with spectroscopic redshitts (z < 4.0) p ublished in Wirth 
respectively. Objects in ICohen et all feOOOT) which are listed in the catalog by I Wirth et al.l [ j2004h 



Wirth et all ll2004h and lCohen et al.l l|2000h 



are omitted. Colour track lines and 



star symbols are the same as in Fig. [2] Filled green circles indicate colours of galaxies at z > 4.5 confirmed by spectroscopy. Three 
open green circles show colours of objects which are not included in our LBG candidates catalog but are reported to be at z > 4.5 in 
previous spectroscopic observations (see table [3J. (b): same as (a) but fo r the J0053+ 1 234 re gion. Squares in magenta are objects with 
spectroscopic redshifts smaller than 4.0 in the reg i on, m ostly taken from lCohen et all jl999aT l. and green circles are galaxies at z > 4.2 
identified by spectroscopy by us and Steidcl et al. (1999) (see table [4}. See text for more details of these figures. 



latter criterion is slightly redder than ours (in both llO j and 
the present paper we use V — I c > 7.0(7 C — z') + 0.15 when 
on the AB normalization), contrary to their argument. In- 
stead, the main cause of the small number counts with their 
selection is that their first criterion sets too strong a lower 
limit in V — I c colour and thus misses many z ~ 4.5 galaxies. 
As is shown in Fig. [3] their criterion, V — 7 C 2.4, misses 
many LBGs with high-z spectroscopic identifications. 

3.6 UV colour distribution 

If our LBG candidates are real objects at z ~ 5 then their 
I c — z' colours represent rest-frame UV colours. In Fig. [5] 
we show the distribution of I c — z colours along apparent 
z'-band magnitude. In this figure we also show the distri- 
bution of I c — z' colours for objects with 23.0 < z' < 26.0 
as filled circles and the median colour values in different z'- 
band magnitudes (in a 0.5 mag step) as squares. Typical 
photometric errors in I c — z' colours are calculated from the 
simulations using artificial objects (described in section 4.1), 
assuming I c — z' = 0.07 mag (from a model SED at z — 4.8, 



E(B — V) = 0.2). Because upper limits on 7 C — z colour val- 
ues are imposed during the selection of LBGs, there are few 
objects with 7 C — z' > 0.4. This does not mean that all star- 
forming galaxies at z ~ 5 have such blue UV colours. We do 
not know what fraction of star- forming galaxies with red UV 
colours (probably due to heavy attenuation by dust) at the 
target redshift are missed in our search; at least in the red- 
shift range z = 2-3 significant populations of star-forming 
galaxies w ith redder rest-frame UV co l our have been de - 
tected (e.g. lForster Schreiber et~al j2004lDaddT et al. Il2004l ). 
Application of similar methods at higher redshifts will be an 
important subject of future space-based infrared missions, 
but for now we must press on with UV-selected samples. 

In the median values of I c — z' colours there is a weak 
trend in that objects bright in z'-band have redder I c — z' 
colours. In other words, the ratio of the number of blue (e.g., 
I c — z < 0) objects to the number of red objects is smaller 
for bright objects than for fainter ones. Since the upper limit 
on I c — z' colour in our selection criteria depends on V — 
7 C colour, the selection effect is not straightforward in this 
figure. However, since there is no restriction for the selection 
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Figure 4. Sky distribution of LBG candidates in the HDF-N region (a) and in the J0053+1234 region (b). Filled magenta circles show 
the positions of LBG candidates. Symbol size indicates objects' brightness in a 0.5 mag step. Filled green circles are for objects which 
have been confirmed to be at z > 4.5 with optical spectroscopy. Large circles with dashed lines indicate the positions and areas of optical 
spectroscopy by us. Gray circles indicate the position of masks used to avoid the effects of bright stars. Rough boundaries of image areas 
used for selection of LBG candidates are shown with dashed lines. Note that not all of the area within these lines is used: due to the gaps 
between CCD chips, different position angles and dither scale lengths, some positions near boundaries have smaller effective exposure 
time, and such areas were not used for candidate selection. Large and small areas drawn with solid lines in Fig. |4f a) are the HST survey 
areas of the GOODS-N and the HDF-N, respectively. 



N 




Figure 5. .z'-band magnitude against I c — z' colour for LBG 
candidates at z ~ 5 in the HDF-N region. Open circles show 
magnitudes and colours of individual galaxies. Typical errors in 
ic — z' colours are displayed in the right side of the panel. Solid 
squares show the median of I c — z' colours in a 0.5 mag step 
in z'-band magnitude. Filled circles indicate the distribution of 
fraction of I c — z' colours in an 0.1 dex step for objects with 
23.0 < z' < 26.0, and its scale is shown in the right side vertical 
axis. 



of objects with blue I c — z colours, the absence of bright 
objects with blue rest-frame UV colour must be real. There 
is no selection bias due to the limiting magnitude in 7 c -band, 
since this is I c — 26.86 mag in the HDF-N region. The red 
colour of UV continua can be attributed to various origins, 
including heavier dust attenuation, older stellar populations, 
and weaker Lyman a emission (for objects at z > 4.8 the I c 
filter covers the rest-frame Lyman a wavelength). We will 
briefly discuss this dependence of the UV colour distribution 
on UV luminosity in section 5.2.2. 



4 UV LUMINOSITY FUNCTION AT Z ~ 5 

4.1 Correction for Incompleteness and 
Contamination 

The procedure that we use to derive the UV lumi nosit y func- 
tion follows that described in our previous paper (|l03t ) which 
was based on shallower data for the HDF-N region only. To 
estimate the volumes sampled with our adopted colour cri- 
teria, we execute a large set of tests in which we use artificial 
galaxies which mimic the size distribution of real galaxies in 
our sample. We are able to derive the completeness of our 
survey for objects with specific magnitudes and colours by 
counting the number of artificial objects which have been 
recovered by the object detection and colour selection pro- 
cedures. In a single test run we generate 25,000 artificial 
objects for each bin of 0.5 in magnitude, 0.1 in redshift, 
and reddening by dust (five different degree of reddening 
from E(B - V) = 0.0 to 0.4 with steps of 0.1). Colours are 
assigned based on the model spectra of star- forming galax- 
ies ge nerated with PEGASE.2 (see Fig. |2), ICalzetti et al.l 
(2000) extinction law, and a prescription for IGM attenua- 
tion bv llnoue et al.l l|2005t h These objects are then put into 
random positions in the V, I c and z' images. The meth- 
ods of object detection and sample selection according to 
the colour criteria are the same as those for real objects, 
described in sections 3.3 and 3.4. 

In Fig. [6] we show the completeness of our survey as a 
function of redshift and z'-band magnitude. Results for five 
different amounts of reddening were averaged with weights 
estimated from the UV colour distribution of real galaxies 
(section 3.6). The maximum completeness is achieved for 
brightest (23.0-23.5 mag) objects at redshift 4.5-5.0 in the 
HDF-N region. We note that there is a significant (> 10%) 
detection rate for objects at z ~ 4. It suggests that in 
our sample of z ~ 5 LBG candidates there should be a 
small number of objects at such lower redshift ranges. So 
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Figure 6. Detection rates of the galaxies for the HDF-N re- 
gion (solid lines) and the J0053+1234 region (dashed lines), as a 
function of z'-band magnitude and redshift. The lines represent 
2'-band magnitude of 23.0-23.5, 24.0-24.5, 25.0-25.5 and 26.0- 
26.5 (the last bin for the HDF-N region only). 

far we have detected one object at z = 4.26 in the optical 
spect roscopy of the J0053+1234 region (tablel4l lAndo et al.l 
120071 ). The existence of such relatively low redshift objects 
does not have a serious effect in deriving the UVLF at z ~ 5, 
since we properly treat them by calculating the survey vol- 
ume using these selection functions. The detection rates in 
the J0053+1234 region are lower than those in the HDF-N 
region at the same magnitude ran ges, and they are similar 
to those for the data used in Il03l . which were ~ 0.8 mag 
shallower in z'-band than the present HDF-N data. 

A weighted average value of redshift derived from this 
completeness distribution is 4.8 for both the HDF-N and the 
J0053+1234 regions. Hereafter we use this value as a typical 
redshift in the presentation of the UV luminosity function. 
From these completeness estimates for different magnitudes, 
redshifts and am ount of dust redden ing, we calculated the 
effective volumes (|Steidel et alj[l999h of our surveys, which 
are listed in Tables [5] and [6] 

As we did in llOSt we use a resampling method to esti- 
mate the rate of contamination by galactic stars and galax- 
ies at intermediate redshifts. We perturbed the colours of 
detected objects which lie outside of the colour selection 
criteria by adding random errors based on the error mea- 
surements in photometry. Then we counted the number of 
objects which fell into the colour selection criteria after the 
addition of these random errors. For each bin of 0.5 magni- 
tude these resampling tests were executed 1,000 times, and 
average numbers of interlopers were derived. These numbers 
are listed in Table and E] as iVjnt-0 

We compute the UVLF as 

*(m) = 2 x (N{m) - AU(m))/Kff(m), 

where N(m) represents the number of LBG candidates de- 
tected within the magnitude bin of m ± 0.25 mag, and 

1 We also tested a case in which we perturbed the colours of 
objects within the colour selection area and count the number 
of objects which get out from the area, and confirmed that the 
changes in the UVLF are smaller t han or comparable to the size 
of overall errors. See section 4.1 of Il03l for more deitals. 



Veff(m) is the effective volume. The multiplication by 2 is 
required to obtain number densities per magnitude, since we 
take 0.5 mag bins in our sample. The UVLFs derived for the 
HDF-N region and the J0053+1234 region are listed in Ta- 
bles [5] and[6] respectively, and they are shown in Fig.[7]with 
open circles (HDF-N) and open triangles (J0053+1234). The 
error bars indicate la error estimates including statistical 
(Poisson) noise, photometric errors and uncertainties in con- 
taminations and completeness. In converting from apparent 
z'-band magnitude to absolute magnitude we assume a lu- 
minosity distance of an object at z — 4.8 under the adopted 
cosmology (Mw = tn z i — 46.33, !Jm = 0.3, Qa = 0.7, 
Ho — 70 km/s/Mpc). In Il03l we made a correction for lumi- 
nosity in the conversion from 7 c -band apparent magnitude 
to rest-frame luminosity at 1700A, assuming the SED of a 
model star-forming galaxy. In the present analysis we do not 
make such a correction in the conversion from magnitudes 
in z'-band (the effective wavelength is 9086A), and thus the 
LFs we present are at ~ 1570A. 

In the calculation of selection functions and effective 
volumes we assume a dust attenuation distribution from 
E(B—V) = 0.0 to 0.4 based on the I c — z' colour distribution. 
In order to examine the robustness of the LF against this 
assumption about dust attenuation, we also calculated ef- 
fective volumes with two extreme cases in dust attenuation: 
one case is that using just a model without dust attenua- 
tion (E(B — V) = 0.0) and the other is only using a model 
with relatively large dust attenuation (E(B — V) = 0.4). 
In these extreme cases the shape of selection function along 
redshift changes from that shown in Fig. [5] because the red- 
shift range covered with our selection criteria changes when 
the assumed model is different. However, we find that the 
differences in the resulting LF — even with these extreme 
cases — are small, as shown in Fig. [8] In this figure the 
UVLFs in the HDF-N region using effective volumes cal- 
culated based on these extreme assumptions are compared 
with the UVLF computed using the fiducial dust attenu- 
ation distribution. Since the extreme assumptions in these 
tests (all objects are free from dust or all objects have fairly 
large dust attenuation of E(B — V) = 0.4) are unlikely, we 
can say that our UVLF is robust against the assumption of 
model colours and distribution of amount of dust attenua- 
tion. 

The modification of the selection criteria should also 
alter the effective volume. We examined the change of ef- 
fective volumes and luminosity function by means of a test 
in which the colour selection criteria were shifted in both 
V — I c and I c — z with 0.05-0.15 magnitude, and derived 
luminosity functions in the same way as we did for the orig- 
inal selection criteria. The UVLFs obtained by these tests 
show differences smaller than the error bars of our LF shown 
in Fig. [7] With these tests we verified that our UVLF results 
are robust against slight changes to the colour selection cri- 
teria. 

4.2 Combination of the UVLFs in the two survey 
fields and Schechter function fitting 

Since the J0053+1234 sample is shallower than the HDF-N 
sample, the UVLF of the J0053+1234 region is limited to 
z' < 25.5 mag. As shown in Fig. [7] the number densities 
of these two fields agree fairly well, although in two bright 
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Figure 7. UV luminosity functions (UVLFs) at z ~ 5. Open circles and open triangles represent UVLFs of LBGs at z ~ 5 in the 
HDF-N region and the J0053+1234 region, respectively. Filled circles show the UVLF derived by combining results of the two fields. The 
data points of UVLFs based on single fields are shifted +0.05 mag for clarity. The solid line represents the Schechter function fitted to 
the combined UVLF, and the dashed line is a fitting for the UVLF in the HDF-N. Open squares are the UVLF in the HDF-N region 
determined with previous shallower data (103). 
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Figure 8. The UVLF at z ~ 5 for the HDF-N region (same as the UVLF with open circles in Fig. calculated using the distribution 
of dust attenuation estimated from I c — z' colours (Fig. [5} is compared with UVLFs derived based on two extreme assumptions about 
the dust attenuation of our sample LBGs (no attenuation and E(B — V) = 0.4 for all objects), shown with open circles. LFs with open 
circles are shifted +0.05 mag for clarity. See section 4.1 for details. 
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magnitude bins (23.5 < z' < 24.0 and 24.0 < z' < 24.5) 
the discrepancy between results from the two fields is large 
enough that it exceeds the size of the 1-a error bars. This is 
probably due to the small number of objects (5-35) in these 
magnitude bins, and it might also indicate the existence of 
field-to-field variations, the effect of which is not included in 
the error estimates. 



We averaged the number densities of the two fields by 
applying weights based on their survey areas (for magnitude 
range brighter than z — 25.5), and derived the final UVLF 
with the present data set. This UVLF is listed in Table [7] 
and is shown as filled circles in Fig. [7] Hereafter we call this 
the final UVLF. 



The results of fittin g the UVLFs with the Schechter 
function jSchechtedl 19761 ) are also shown in Fig. [7] The solid 
line is a function fitted to the final UVLF, and the dashed 
line is for the UVLF based solely on the deeper HDF-N 
region data. The best-fitted parameter values for the final 
UVLF and for the LF from the HDF-N region are summa- 
rized in Table E] We note that the z ~ 5 Schec hter param- 
eter fit given in Saw icki and Thompson! l|2006bD is based on 
a (slightly earlier) version of the HDF-N region UVLF pre- 
sented here. Fig. also shows the 68% and 95% confidence 
level error contours calculated via Monte-Carlo resampling 
of the UVLFs. For each point of the observed UVLFs we 
add a random error according to its uncertainty (assuming 
Gaussian distribution), and execute the Schechter function 
fitting to this perturbed UVLF and obtain a \ 2 value. We 
repeat this test 1,000 times and derive the region in param- 
eter space that contains 68% or 95% of these resampling 
results. 



Although the number densities in the combined data 
and those in the HDF-N region are fairly close in the fainter 
part of the LF (24.0 < z < 25.5), the best-fitted Schechter 
function parameters are significantly different. This differ- 
ence is caused by relatively large difference in the number 
density in the bright 23.5 < z' < 24.0 magnitude range. 
Even the faint-end slope a is affected by the change in the 
LF's bright part, due to changes in <f>* and M* . The large 
field-to-field variation in the bright end of the UVLF be- 
tween our two survey areas (both fields have areas larger 
than 500 arcmin 2 ) clearly indicates the necessity of a large 
survey area and multiple survey fields. A study based on a 
single field observation may be affected by large-scale struc- 
ture in the galaxy distribution. Moreover, if the survey area 
is small, it is very likely that the survey fails to detect bright 
and rare objects and so underestimates their number den- 
sity. We also note that the Schechter function parameters 
are not a good way to represent the properties of observed 
LFs unless an observation reaches far deeper than M* . A 
comparison of LFs at different epochs or in different envi- 
ronments should be cautious of uncertainties and degener- 
acy between parameters (as illustrated in Fig. |5J). The di- 
rect comparison of number densities at in a fixed luminosity 
range is a better way to compare LFs. 



4.3 Comparisons with UVLFs in previous studies 
of galaxies at z ~ 5 

4-3.1 Comparison with the UVLF in Iwata et al. (2003) 

The new UVLFs at z ~ 5 agree well with our previous result 
(|I03T) based on shallower data in the HDF-N region, which 
is indi cated as open squares in Fig. [7] Note that because 
in Il03l we used the Vega-based magnitude system and the 
fiducial redshift was z = 5.0, the magnitude bins are slightly 
different from our present results. 



4.3.2 Comparison with the UVLF of FORS Deep Field 

In Fig. [10] we compare our UVLF with the one at 1500 A 
which is base d on the photometric re dshift sample selected 
in /-band by iGabasch et al.l l|2004ah . These authors used 
150 galaxies in the FORS Deep Field (FDF) with /-band 
magnitudes brighter than 26.8 mag which have estimated 
redshifts between 4.01 and 5.01. While the depth of their 
sample is comparable to ourfl the number of objects is 
small, mainly due to their small survey field area (7' x 7'). 
The ov erall shape of the UVLF derived by IGabasch et al.l 
|2004al ) is similar to ours. The number density of faint galax- 
ies in our UV LF at z ~ 5 is s lightly smaller, as has been 
pointed out by G abasch et al.l (|2004aT ) who used our previ- 
ous UVLF for the compa r ison. S ince the redshift range of 
galaxies in IGabasch et all (|2004al ) is smaller than ours, this 
trend might indicate a real number density evolution of UV- 
selected star- forming galaxies (see section 5.1 for discussion 
of the differential UVLF evolution), although the difference 
in the methods of sample selection and possible field-to-field 
variation prevent us to make further discussions. 



4-3.3 Comparison with the result of Subaru Deep Survey 

In Fig. QJj] we show with crosses t he UVLF of LBGs at z ~ 
4.7 derived bv lOuchi et alj (1200481 ) in the Subaru Deep Field 
and show their fitted Schechter function using the dotted 
line. Note th at their faint - end slo pe a was fixed to —1.6. As 
described in lOuchi et al.1 (l2004al ). there is a discrepancy in 
the bri ght p art (Muv < —21.5 mag) between our previous 
UVLF |l03h and their result. The discrepancy still remains 
when o ur updated data is compared with the lOuchi et al.l 
(|2004ah UVLF, and this seems somewhat unnatural, since 
the field-to-field variance between our two survey fields (at 
M\jy < —22) is more than three times smal ler than the 
differen ce betwee n our UVLFs and the one bv lOuchi et al.l 
|2004ah . Recently. lYoshida et al.1 12006) derived an updated 
UVLF of z ~ 5 LBGs in the same Subaru Deep Field using a 
deeper and «160 arcmin 2 wider image data set. Their filter 
set, procedures for reducing data, selecting LBG candidates, 
and calculat i ng the UVLF are basically the same as those in 
lOuchi et al . (2004a). Conseq uently, their results agree with 
those bv lOuchi et al.l l|2004ah and thus the number densities 
in their LF's bright part remains significantly smaller than 



2 We should note th at the g- and R- band images used in 
IGabasch et al. (2001a) have depths comparable to that in /-band 
and may not be deep enough to detect the spectral break reliably, 
while their B-band image is ~ 0.8 mag deeper than /-band. 
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Figure 9. Confidence contours of Schechter function parameters. Solid lines show areas show the 68% (inner) and 95% (outer) confidence 
levels for the 'final' UVLF which combines the results of the two fields. Dashed lines are the same as the solid lines, but for the (deeper) 
HDF-N data only. Small points show the values of the best-fitted parameters. 



2 AB mag (for 2 = 4.8) 
23 24 25 26 27 

- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 |^ .1 - ■ — - H — - I 




-23 -22 -21 -20 -19 

M AB (UV) 

Figure 10. UV luminosity functions at z ~ 5. Filled circles show our UVLF derived by combining results of the two fields. This is 
the same UVLF as in Fig. The solid line repres ents the Schech t er func tion fitted to this combined UVLF. Solid squares represent the 
rest-frame 1500A luminosity function derived by iGabasch et al.l ll2004ar> from their /-band selected photometric redshift sample with 
estimated redshifts at 4.01 < z < 5.01, and the dot-d ashed line is a Schechter function fitted to it. The crosses are the UVLF for LBGs 
at z ~ 5 (U-band drop LBGs) bv lOuchi et al. | ll2004ah and the dotted line represents the Schechter function fitted to their data with the 
faint-end slope fixed to a = —1.6. 
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We further investigated the origin of the differen ce be- 
tween our results and those by lOuchi et al.l |2004a). One 
possible cause of this discrepancy is the difference in colour 
selection criteria. However, we must ta ke care here because 
the filters use d by us and those us ed bv lOuchi et alj (|2004al ) 
are different: lOuchi et al.l (|2004al ) used i' while we used 7 C , 
and since the i'-band filter has transmission that is shifted 
slightly to a shorter wavelength with respect to the J c -band 
(see Fig. [TJ, the comparison of sample selection criteria is 
not straightforward. We illustrate this problem as follows. 
In Fig. [11] we use a dot-dashe d line to show the V — i' and 
i' — z colour criteria used in lOuchi et all (|2004al ). Fig. [TT] 
also shows the redshift tracks of star-forming galaxies in the 
colour pair V — I c and I c — z' (solid lines) and in the colour 
pair V — i' and i' — z' (dashed lines). Since the effective wave- 
length of i'-band is slightly shorter than that of 7 c -band, 
the i' — z colour becomes red at relatively low redshif t (at 
z ~ 4.7). If we were to apply the lOuchi et all (|2004al ) cri- 
teria directly to our V — I c and I c — z' colours, then the 
sample would be contaminated by many low- z galaxies and 
galactic stars with V — I c > 2.4 and I c — z' > 0.4, as shown 
with red squares and star symb ols in Fig. | llfa). Ho wever, 
the selection criteria adopted bv lOuchi et al. (|2004al ) avoid 
including these objects because the i' — z' colours of these 
objects are ~ 0.2 mag redder than their I c — z' colours. This 
illustrates the difficulty in making a direct comparison of 
the colour selection criteria of the two surveys. 

Thus, in order to roughly estimate the effect of the dif- 
ference between colour selection criteri a, we adopt modified 
criteria which consist of those used bv lOuchi et al.l (|2004al ) 
at I c — z' < 0.3 and those used by us at I c — z ^ 0.3. 
These modified criteria are shown with a thick solid line in 
Fig. Illf b); they approximately reproduce the lOuchi et al.l 
(20041) criteri; as applied to our V — I c , I c — z' data. 
In Fig. [12] we show z'-band number counts of LBGs in 
our survey field s, using bot h our c riteria and the mod- 
ified criteria of lOuchi et al.l ll2004al) and compare them 
with the number co unts given by lOuchi et al.l (|2004al ) and 
lYoshida et al. (2006). The change of colour selection criteria 
reduces the number of LBGs in the HDF-N region, but there 
remains a clear discrepancy betwee n our number counts and 
those in the Subar u Deep Field l|Quchi et al.l l2004al and 
lYoshida et ail [2006). In the 875 arcmin of Yoshi da et al.l 
1 20061 ) there is only one object with z' < 24.05 mag and only 
12 objects with z:'-band magnitude between 24.05 and 24.55, 
while there are 7 objects with z < 24.0 mag and 35 objects 
with 24.0 < z < 24.5 mag in the HDF-N region when we 
use our own selection criteria. As shown in Fig. Illf b) , when 
we adopt the modified colour criteria, six out of 7 LBG can- 
didates with 23.0 < z < 24.0 remain within the selection 
area. In 24.0 < z < 24.5 the number of objects becomes 
17, about half of our original sample, but th i s corr esponds 
to ~ 29 for the survey area of lYoshida et alj (|2006T l and so 
represents a surface number density more than two times 
larger than theirs. Thus, the modification in the colour se- 
lection criteria is unable to explain the differences between 
the number densities in our final UVLF and those in the 
Subaru Deep Field, which are 2-5 times larger than the la 
errors estimated for our UVLF. 

In addition to the change in the number density of ob- 
jects, the change of colour selection criteria reduces the ef- 
fective volume. These two effects counteract each other and, 



as a result, the luminosity function would remain relatively 
unchanged. Because of the difference in the filter sets, we 
cannot calculate the effective volu me in the cas e of app ly- 
ing the colour selection criteria of lOuchi et alj (|2004al ) to 
our sample. We are thus unable to calculate the UVLF with 
such criteria and so must remain focused on comparing raw 
number counts. However, as we have mentioned in section 
4.1, slight changes to the colour selection criteria do not 
change the shape of the UVLF. 

Although the number of objects in our sample which 
were observed with optical spectroscopy is still limited, 6 
out of 8 objects with z'-band magnitude between 24.0 and 
24.5 mag in the HDF-N region for which we obtained spec- 
troscopy are confi rmed to be at z > 4.5 (including one AGN; 
lAndo et alj|2004 ). The redshifts of the remaining two ob- 
jects are unidentified due to the absence of clear spectral 
features. Thus, we believe that most of bright objects in our 
sample are genuine LBGs at z ~ 5. Indeed, if the success 
rate of spectroscopic confirmation (6/8) for bright LBGs is 
representative of that for the whole of our sample, it im- 
plies that at least 26 among the 35 HDF-N-region galaxies 
with z' = 24.0-24.5 are at z > 4.5. This number density of 
26/505.85=0.051 arcmin -2 well exceeds the number density 
of candidates in the same mag n itude range (12/875=0.013 



N /arcmin 2 ) in Yoshida et al. (|2006l). Consequentl y ; it is 



possible that lOuchi et al. ( 20043) and lYoshida et al] (j2006h 
may miss significant numbers of luminous high-z objects un- 
less the field-to-field variance between the HDF-N region 
and SDF is far larger than expected. However, such large 
field-to-field variance seems unrealistic considering the spa- 
tial distribution of LBG candidates in our survey fields and 
the similarity of the UVLF results in our two fields. 

Thus, although the cause of this discrepancy remains 
an open question, because of the comparisons and consider- 
ations described above, we think our UVLF is robust. Here- 
after we use our final LF result as a representative z ~ 5 
UVLF to compare with UVLFs at lower redshift ranges. 



4.4 Comparisons with results of numerical 
simulations 

Recent development of computer resources and numerical 
techniques makes comparisons of observational results with 
results from numerical simulations interesting since such 
comparisons make it possible to investig ate the phys i cs be - 
hind the observed properties of galaxies. iNight et all (|2006l ) 
executed cosmological SPH simulations, calculated SEDs for 
galaxies in their simulations using the population synthe- 
sis code of iBruzual and Chariot! (|2003l ) , and then applied 
LBG colour selection criteria to create a simulated sample of 
LBGs at z =4-6. They use this sample to compare UVLFs in 
their numerical si mulations with ours ( UVLF for the HDF- 
N regio n) and withlOuchi et alj |2CI04a) . Since the prescrip- 
tions of INight et alj l|2006l ) do not include attenuation by 
dust, their UVLFs can be shifted along the luminosity axis 
by varying the assumed amount of uniform dust attenua- 
tion. The UVLF of z = 5 LBGs in their simulations can be 
m ade to agree with o bserved UVLFs (both ours and that 
bv lOuchi et ail2004al ) by varying dust attenuation between 
E(B — V) = 0.15 a nd 0.3 (using t he at tenuation law by 
ICalzetti et al. 2000). lFinlator et all (|2006l ) applied dust at- 
tenuation based on a correlation between the metallicity and 
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Figure 11. V — I c and I c — z' two colour diagrams, (a): Red and blue squares are galaxies at z < 4.0 in the HDF-N region, same as in 
Fig. [3^.. Green circles are spectroscopically confirmed objects at z > 4.5 in the HDF-N. For objects not detected in F-band, arrows are 
used to show the 3-cr lower limits of V — I c colours. The dashed line and shaded area stands for colour selection criteria adopted in this 
study for the selection of Lyman break galaxies at z ~ 5. The colour selection criteria of Ouchi et al. (2004a) (used with V — i' and i' — z' 
colours) are displayed with a dot-dashed line. Blue lines are colour tracks of a star-forming galaxy without dust attenuation, and lines 
in orange are for dust attenuation of E(B — V) = 0.4. Solid lines show colours in V — I c and I c — z' , and dashed lines are in V — i' and 
i' — z' . Star symbols indicate the V — I c and I c — z' colours of A0-M9 stars calculated based on the library by Pickles (1998). (b): same 
as (a), but shows objects with 23.0 < z' < 24.5 mag detected in our V, I c and .z'-band images for the HDF-N region. Sizes of symbols 
indicate object brightness in 0.5 mag bins, with the largest symbols representing objects with 23.0 < z' < 23.5 mag. The thick solid line 
represents the modified colour selection criteria used for the test described in section 4.3.3., which consist of the colour selection criteria 
of Ouchi et al. (2004a) for i' — z' < 0.3, and for »' — *'> 0.3 follows our colour selection criterion. 



dust attenuation inferred from the large sample of nearby 
galaxies taken from the Sloan Digital Sky Survey. They 
only show galaxies at z = 4 in their numerical simulations, 
and their UVLF is broad ly consistent with obse r vationa l 
results at z ~ 4 su ch as Sawicki and Thompson! (|2006al ). 
lOuchi et al] l|2004al ) and iGabasch et aD (l2004a(). Although 



the U VLFs found in num e rical s imulations by lNight et al] 
(2006) and iFinlator et all (|2006l ) show fairly steep slopes, 
a > —1.7, such slopes are estimated from magnitude ranges 
fainter than the limiting magnitudes of current deep obser- 
vations (m > 27). Thus, in summary, no serious conflict 
with observed results has been found in these numerical 
simulations. However, we should note that there are free 
parameters in models used in these simulations which are 
not strongly constrained by physics in a self-consistent way 
(e.g., amount of dust attenuation, effects of feedback from 
star formation and AGN), and it is still not possible to dis- 
criminate the "best" observed UVLF through the compar- 
isons with results of such numerical simulations. One inter- 



esting prediction common in these two simulations is that 
UV-luminous objects should have higher stellar mass than 
low-luminosity ones. This trend may be related to the dif- 
ferential evolution of galaxies we describe in section 5, and 
this prediction is testable with deep infrared observations of 
LBGs at z > 4. 



5 DIFFERENTIAL EVOLUTION OF UV 

LUMINOSITY FUNCTION FROM Z ~ 5 TO 3 

5.1 Comparison with the z ~ 4 and 3 UVLFs by 
Sawicki and Thompson (2006) 

Sawicki and Thompson i|2005l ; l2006al ) constructed ultra- 
deep samples of UV-selected star-forming galaxies at 1.7 < 
z < 4 using the very same U n G1ZI filters and colour selec- 
tion criteria as those used at brighter magnitudes by Stei- 
del and his collaborato r s. Usin g these very deep samples, 
ISawicki and Thompson! l|2006al ) derived the UVLFs reach- 
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Figure 12. The z'-band number count of LBGs at z ~ 5. Filled circles and triangles are those for the HDF-N and J0053+1234 regions, 
respectively, based o n the sample selecte d using our colour cr iteria. Number counts of V-drop LBGs in SDF are shown with open circles 
llOuchi et al.ll2004al) and open triangles dYoshida et al.ll2006f) . Filled squares show the number counts of LBGs in the HDF-N selected 
using the modified colour criteria based on those used in Ouchi ct al.(2004). See text for details. No correction for incompleteness has 
been made. Vertical error bars show simple Poisson errors. 



ing far below L* z=3 . Specifically, at z — 3 their limiting mag- 
nitude is Mu m — —18.5 mag and Mu m is about —19.0 mag 
at z — 4. These depths are about 1.5 magnitudes fainte r 
than previous large-area studies (e.g., St eidel et al] fl999). 
but the results are highly reliable because of the use of the 
selection technique identical (down to the same filter set 
and color-color selection criteria) to that in the Steidel et 
al. work. Comparisons of our new z — 5 UVLFs with these 
deep UVLFs at lower redshifts can lead to crucial insights 
regarding the evolution of UV-selected star-forming galaxies 
during early epochs in the history of galaxy formation. 

In Fig. [13] we show our new z — 5 UVLF (the final 
weight-averaged LF, as described in sec. 4.2) along with 
the z — 4 and z — 3 UVLFs from ISawicki and Thompson! 
(2006a|). In ll03l we compared the previous, shal lower version 
of our 2-5 UVLFs with z = 3-4 UVLFs bv lSteidel et~al] 
i| 19991 ) and argued that there is no significant difference be- 
tween UVLFs from z = 3 to 5, although there might be a 
slight decrease in the fainter part. This suggested trend now 
becomes clear thanks to the availability of deeper galaxy 
samples at all three epochs: although the limiting absolute 
magnitude of our z — 5 sample is shallower than the lower- 2 
samples of ISawicki and Thompson! (|2006al ). we see that at 
luminosities fainter than L* =3 (Muv ~ —21.0) the number 
density of LBGs steadily increases from z — 5 to 3, while no 
number density evolution is found at brighter magnitudes. 

In order to be rigorous we need to quantify the strength 
and statistical significance of this differential, luminosity- 
dependent evolution. However, we do not wish to compare 
the evolution of the Schechter function parameters because 
the three Schechter parameters (M* , (f>* , and a) are highly 
degenerate and so it is difficult to physically interpret the 



meaning of their time-evolution. Instead, it would be far 
better to study the evolution of the observed number den- 
sities directly, but this is not simple since, at the three dif- 
ferent epochs, the UVLF data span somewhat different lu- 
minosity ranges, and because the magnitude bins are cen- 
tred at different Mtjv- To overcome these problems , we use 
the approach developed bv lSawicki and Thompson] (|2006a) 
— we investigate how far the observed number densities at 
the three epochs deviate from a fiducial model, for which 
we adopt the Schechter function fit to their 2 — 3 data. 
In this way, for each data point, we construct $(Af, z) = 
<f>data{M, z) I '(j)f~ t 3 (M), which measure the departure of the 
data from the model. At z ~ 3, $(M) ~ 1 at all M, re- 
flecting the fact that the 2 — 3 Schechter function is a good 
description of the 2 — 3 data. At other redshifts, however, 
the <I>(M) deviate from unity, reflecting the UVLF's evolu- 
tion. Next, at each epoch we compute $(2) — i.e., the av- 
erage of the <E>(M, 2) — for bins brighter (and fainter) than 
Muv = M*= 3 = —21.0 mag. Finally, the ratio of the <l at 
two redshifts tells us the amount of number density evolu- 
tion that a given population (bright or faint) undergoes. 

The results of the above procedure are shown in Fig. [14] 
(see also Table [SJ, which shows the number density evolu- 
tion of the bright end (horizontal axis) and the faint end 
(vertical axis) of the UVLF. The data show strong evidence 
for evolution of the faint end of the UVLF: at 2 — 5 there 
are only 0.34 ± 0.02 as many faint galaxies as there are at 
2 — 3. This represents a three- fold increase that is highly 
statistically significant. The evolution over 2 — 5 — > 4 is 
smaller (it is a 1/0.78 = 1.3-fold increase), but since the 
cosmic time between 2 = 4.8 and 2 = 4 (—300 Myr) is half 
of the interval between 2 = 4 and 2 = 3 (—600 Myr), the 
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data are completely consistent with the idea that the num- 
ber of faint LBGs increases at a constant rate from z ~ 5 
to 3. However, while the faint end of the UVLF is evolving, 
there is no evidence for any change in the number density of 
luminous (Muv < —21.0 mag) galaxies over z ~ 5 — > 4 — > 3. 
Indeed, the z ~ 5 — * 3 data rule out a scenario — indicated 
by the diagonal line — in which the number density of galax- 
ies evolves uniformly at all luminosities. These z ~5-based 
results are in excellent agreement with the tr ends found at 
z ~ 4 — > 3 bv lSawicki and Thompson! |2006a|) (shown with 
green symbols) , although our new z ~ 5 data both lengthen 
the cosmic time baseline and strengthen the statistical sig- 
nificance of those earlier conclusions. 

In summary, the c omparison of our UVLF at z ~ 5 with 
those at z = 3-4 from lSawicki and Thompson! l|2006al ) indi- 
cates that the evolution of LBG number densities depends 
on UV luminosity. Over the ~1 Gyr from z ~ 5 to z ~ 3, the 
number density of luminous LBGs stays almost unchanged 
while there is a gradual and steady build-up in the number 
density of faint (Muv > —21 mag) galaxies. 



5.2 Implications of UV Luminosity-dependent 
Evolution 

5.2.1 Star Formation History of LBGs 

The UV luminosity-dep endent evolution of UVLF has b een 
already pointed out by ISawicki and Thompson! (l2006al ) for 
their LBG samples at z ~ 3 and 4. Our data indicate that 
similar evolutionary trends were already in place at earlier 

cos mic times. 

ISawicki and Thompso n (2006a) discuss three possible 
scenarios responsible for the differential evolution of the 
UVLF between z ~ 4 and 3. The first one is the straightfor- 
ward increase of the number of low-luminosity LBGs at later 
epochs. They argue that, if the ages of low-luminosity LBGs 
at z ~ 3 are mostly smaller than 600 Myr, the increase in 
number density can be a natural consequence of the emer- 
gence of large population of low-luminosity galaxies after 
z ~ 4. Sufficiently young ages of z ~ 3 LBGs have been sug - 
gested by several authors. For example. IShaplev et al.l (12001 ) 
used rest-frame UV to optical SEDs of their (fairly bright) 
sample LBGs for comparisons with template SEDs based on 
models with a simple constant star formation history and 
derived a median age of 320 Myr, and results with similar 
or younger ag e estimates have been also reported by other 
authors (e.g.. ISawicki and Yed Il998l ; iPapovich et all l200ll ; 
llwata et al.l 120051 ). We should note, however, that the ages 
estimated by SED fitting based on simple star formation his- 
tory models (constant or exponentially decaying SFH) are 
strongly affected by recent star formation activity. Thus the 
age should be interpreted as time since the onset of most re- 
cent star formation, and such age estimates may be incorrect 
if LBGs hav e sp oradic star formation histories. As we have 
discussed in Il03t since the number density of LBGs in the 
luminous part of the UVLF shows little change from z ~ 5 
to 3, the small ages of LBGs at 2 ~ 3 would imply that the 
star formation histories of LBGs are sporadic, and that in- 
dividual LBGs would have experienced short bursts of star 
formation during their evolution. Thus, if the evolution of 
the fainter part of the UVLF represents the emergence of a 
large population of low-luminosity LBGs at later epochs, it 



should be interpreted as the increase in the number density 
of low-luminosity population as a whole, and for individual 
galaxies a low-luminosity LBG at z ~ 4 can be a luminous 
one at z ~ 3 and vice versa. 

The change in the properties of starbursts, such as in- 
tensity, duration and duty cycles, can be the second possi- 
ble cause of differential evolution. It would be possible to 
explain the increase in the number density of faint galax- 
ies by introducing the differential change in the properties 
of starburst depending on U V luminosity (scenario B of 
ISawicki and Thomp son 2006a). If intrinsically less luminous 
galaxies (unexplored with current observations when in their 
low star formation state) progressively spend longer time in 
the relatively luminous phase, they can be observable by 
z ~ 3, and as a consequence the number density in the faint 
end of the UVLF could increase. The luminous end would 
not change if the starburst properties of luminous galaxies 
remain the same. Such scenario might be tested through 
the comparison of starburst ages of LBG populations at 
different epochs. Exploration of the stellar populations of 
LBGs at z > 4 has become possible recently thanks to the 
availability of rest-frame optical data obtained through mid- 
infrared imaging with the Spitzer space telescope, and sev- 
eral authors have reported the existence of massive galax- 
ies at z ~ 6 (e .g., lEvles et al.l 120051 : iMobasher et al.l 120051 : 
lYan et all l2u06^ . Measurement of stellar masses would to 
some extent resolve the ambiguity for the star formation 
history of LBGs present in analyses based solely on the 
rest-frame UV wavelengths. A systematic study based on a 
uniform selection scheme across the redshift range z = 3-5 
would be an interesting project to pursue in order to address 
these issues. 

One of the major obstacles for precise estimation of 
LBGs' ages would be the notorious degeneracy between 
age and dust in their effects on galaxy SEDs. Indeed, 
ISawicki and Thompson! (|2006a! ) raised the change of dust 
properties as the third possible cause of differential evolu- 
tion. If the property of dust changes differentially depend- 
ing on UV luminosity, this would also alter the shape of the 
UV luminosity function. Unfortunately, we know quite little 
about the properties of dust in high-z star- forming galax- 
ies, a nd usually one attenuation law (such as lCalzetti et al.l 
2000) is assumed for all galaxies in a high-z sample. Thus, 
it could be difficult to discriminate whether it is age or 
dust that changes by studying the discrete SEDs obtained 
through the broad-band imaging data. 

5.2.2 Biased Galaxy Evolution 

We recently found that Lyman a equivalent widths of z ~ 5- 
6 LBGs (a nd Lyman a emit ters) tend to depend on UV 
luminosity (lAndo et al.ll2006T ): high luminosity LBGs show 
little or no Lyman a emission, while there are objects with 
Muv > —21 mag which show strong emission lines. Along 
with the result presented in this paper that the number den- 
sity of luminous LBGs at z ~ 5 is almost the same as at later 
epochs, the UV luminosity dependence of Lyman a equiva- 
lent widths is another indication that UV luminous galaxies 
have different pro perties from fainter ones. As discussed in 
lAndo et al.l (|2006T l. there are several possible causes of the 
luminosity dependence of the strength of Lyman a emis- 
sion, one possible explanation being that luminous LBGs 
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Figure 13. UV luminosity function (UVLF) at z ~ 5 (solid circles and solid line) in this study and those at z ~ 4 (dashed line) and at 
2 ~ 3 (dot-dashed line) by Sawicki and Thompson (2006a). 



are embedded in dense reservoirs of neutral hydrogen gas. 
The existence of such large reservoirs of neutral gas could be 
the cause of the large star formation rates we see and also 
could be the underlying cause of the heavy absorption of 
Lyman a emission by ionized gas surrounding the star form- 
ing regions within galaxies. Another possibility is the large 
amount of dust and/or metals in UV- luminous LBGs which 
have experienced significant amounts of star formation until 
z ~ 5. The existe nce of strong silico n and carbon interstellar 
absorption lines (|Ando et al l 2004) in the luminous objects 
in our sample of z ~ 5 LBGs supports this idea, although 
sample size is still quite limited. A third indication of the 
UV luminosity dependence is the relatively small number of 
bright objects with blue rest-frame UV colours in our sample 
(section 3.6 and Fig. [SJ. The dependence on UV luminosity 
of the clustering strength of z = 3-5 LBGs has also been 
reported by several authors (e.g.. iGiavalisco and Dickinson! 
l200ll ; IOuchi et al 1l2004bl ; lKashikawa et alj|2006h . It suggests 
UV luminous LBGs reside in relatively massive dark matter 
haloes. 

The existence of a number density of UV luminous 
LBGs at z ~ 5 that is comparable to those at z = 3-4, 
as well as the luminosity-dependent properties of LBGs es- 
pecially at z ~ 5, support the idea of biased galaxy forma- 
tion. The luminous galaxies, presumably hosted by relatively 
massive dark matter haloes which are also reservoirs of huge 
amounts of neutral gas, might start bursts of star formation 
at early epochs (i.e., z > 5). On the other hand, the num- 
ber density of faint objects (Mtjv > —21 mag) gradually 
increases along with the increase in the number density of 
dark matter haloes. The UV luminous objects, which have 
SFRs as large as 100MQ/yr -1 at z ~ 5, may become fainter 
at later times because the duration of such enormous bursts 
of star formation can be on order of 100 Myr at most given 



that significantly longer bursts would result in masses that 
exceed 10 11 Mq by z ~ 3. However, some fraction of objects 
that are relatively faint at z ~ 5 would be ignited during 
the period between z ~ 5 and 3, and this would make it 
possible to keep the number density of UV-luminous LBGs 
as a whole almost constant over the intervening ~1 Gyr. 
Such mixing of LBG populations that represent different 
generations across cosmic time may by z ~ 3 extinguish the 
dependence of dust / metal or of gas mass on UV luminos- 
ity, and this mixing may explain the apparent absence of 
UV luminosity dependence of Lyman a equivalent widths 
and strengths of interstellar metal absorption lines at these 
lower redshifts l|Shaplev et alj|2003j ). This biased formation 
hypothesis could be tested by further analyses of chemi- 
cal compositions and stellar populations for LBGs at differ- 
ent redshifts and luminosity ranges. Alth ough attempts to 
tackle these issues have been started (e.g. lAndo et al.ll2004j ; 
lEvles et al1l2005l ; ILabbe et all 120071 ). deep follow-up obser- 
vations in both spectroscopy and multi-wavelength imaging 
that cover many sample galaxies are required. 



5.3 Comparison with UVLFs at z > 6 

Recently, there have been several attempts to determine 
the UVLF at z ~ 6 through d e ep o b servations in 
i an d z bands (e.g., Bouwens et al.l 20031; Bunker et al.l 
2004 I Yan and Windhorstl 12004 IShimasaku et alj 120051 ). 
Bouwens et al. made a compilation of very deep ob- 

servations with HST (GOODS-N, GOODS-S, the Hubble 
Ultra Deep Field (UDF) and UDF pararell ACS fields) and 
derived the z ~ 6 UVLF at 1350A down to Mi 350 ~ -17.5 
mag. They argued that there is a strong evolution in the 
number density of luminous (M < —21 mag) LBGs from 
z ~ 6 to z ~ 3, while the number density of faint galax- 
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ies is comparable to or even larger than that at z ~ 3. They 
claim that it could be caused by the brighten i ng of M * along 
cosmic time. Also, iBouwens and Illingeorthl (|2006l ) utilized 
deep HST/NICMOS observations to find LBGs at z « 7-8 
and concluded that there is a deficiency of luminous galax- 
ies at z ss 7-8 compared to z ~ 6, suggesting that the rapid 
evolution of the number density of luminous galaxies has 
started at 2 > 7. Their arguments are opposite from the 
differential evolution of UVLFs we found for the redshift 
range 3-5, and these two sets of results would be inconsis- 
tent if one assumes that the number density evolution from 
z ~ 6 to 3 proceeds continuously. However, we should note 
that the number densities in the fainter part of the z ~ 6 
UVLF reported by several author s have a fairly large scat- 
ter (see Table 13 and Figure 14 in lBouwens et aUl2006l for a 
summary), even though most of these studies are based on 
similar data sets taken from deep HST surveys. Also, cur- 
rent deep HST based studies for z ~ 6 LBGs suffer from 
severe limit ations in survey volum es. The GOODS survey 
area which IBouwens et al.l (|2006l ) relied on to correct for 
the field-to-field variations is just a fourth of the effective 
survey area presented here. Thus the z ~ 6 UVLF and its 
Schechter parameters would be still need to be further ex- 
amined. Although the number of objects with zgso < 25.5 
mag (which roughly corresponds to M£- 3 ) detected in the 



GOODS area by IBouwens et all (|2006h (8 objects) is incon- 
sistent with no evolution from z ~ 5 even considering the 
dimming from z ~ 5 to 6 and the difference in the rest-frame 
wavelength traced by the z-band filter (which may lead to 
~ 0.1 mag of dimming), it is still premature to draw a con- 
clusion that whether the UVLF evolution from z ~ 6 to 3 is 
due to the brightening of L* . Further studies of z ~ 6 LBGs, 
especiall y those based on wide r survey area are indispens- 
able fcf. lshimasaku et alj|2005h . 



6 UV LUMINOSITY DENSITY 

6.1 UV Luminosity Density at z ~ 5 

Using the updated luminosity function of z ~ 5 LBGs we 
can calculate the UV lu min osity density at that epoch more 
reliably than we did in Il03l thanks to the deeper and wider 
data. The luminosity density at a target epoch due to light 
emitted from galaxies brighter than Li im can be calculated 
by the integration of the LF as 



p(L > Llin 



L</>(L)dL. 



We take two ranges of integration. One lower limit in the 
integration is Mtjv ^ —18.5, which corresponds to 0.1L*_ 3 , 
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and the other is Lu m = 0. The formal lower lu minosity limit 
has been used by several previous studies (e.g., Steidel et alj 
1 19991 : ISawicki and Thompson! l2006bl : iBunker et al.l l2006f ). 
and correspond s to the limiting magni t ude of the 2 ~ 3 
LBG sample of ISawicki and Thompsonl l|2006al ). At z ~ 5, 
we have to extrapolate our observed UVLF more than 1 
magnitude in order to obtain the UV luminosity density 
integrated down to Afuv ^ —18.5 mag. With the extrap- 
olation of the integral down to zero luminosity [Lu m = 0) 
we want to estimate the total UV luminosity density emit- 
ted from star-formation galaxies into inter-galactic space. 
By comparing the UV luminosity densities down to 0.1Ll =3 
with the total ones we are able to estimate the contribution 
from faint objects. The steepness of the UVLF faint-end 
slope has a crucial importance for the total UV luminosity. 
As we have stated in section 4.2, the faint-end slope of the 
UVLFs is still not well constrained by observations (espe- 
cially for higher redshift, i.e., z > 4), and it introduces a 
large uncertainty in the estimation of total UV luminosity 
density. 

Fig. fK] shows with filled circles the UV luminosity den- 
sities obtained from our data at z = 4.8. The integrations 
are made using the best-fit Schechter function for the final 
UVLF described in section 4.2, and error bars indicate la er- 
rors calculated using Schechter functions for perturbed data 
generated to obtain 68% confidence parameter ranges (see 
section 4.2). Table ITTJ1 lists the UV luminosity densities at 
2 ~ 5 calculated using the final UVLF and the UVLF for the 
HDF-N region. Our survey only reaches to Muv < —19.8, 
and thus 68% (49%) of UV luminosity density integrated 
down to Muv ^ —18.5 (Lu-m = 0) comes from this observed 
luminosity range. Since our upd ated UVLF presented here 
is consistent with our previous (|l03l) one in the luminosity 
range covered by the previous data, the UV luminosity den- 
sity at z ~ 5 also shows little change from the value in Il03l . 

6.2 Evolution of UV Luminosity Density 

In Fig.[T2]we also show results for various redshifts by other 
authors. For z ^ 1, we show UV luminosity density data 
from the co mbined sample of G ALEX and VIMOS-VLT 
Deep survey ([Ta kcuchi ct Iil1l2005l . who took the UVLF data 
from lArnouts et al.l 20051 ) . The UV luminosity de n sities a t 
z ~ 2, 3 and 4 are from ISawicki and Thompsonl (2006bl). 
Their r esults are based on UVLFs in lSawicki and Thompsonl 
(1200631 ). Redshift 4.0 and 4.7 results bv lOuchi et~ai1 (1200481 ) 
are also shown. At redshift 6 w e show results by two dif- 
ferent authors: triangl es are from Bouwcns et al. I |2006l ) and 
the pentagon is from IBunker et al.l (|2006l ) (only the value 
for L 0.1L*= 3 is provided in the latter paper). The rest- 
frame wavelengths of all of these observations lie within 
1350A-1700A. For these published values, we only show er- 
ror bars when errors for these specific integration ranges 
(L ^ 0.1L* =3 and L > 0) are provided in the papers, since 
we do not know the range of Schechter LF parameters for 
la confidence level which is required to precisely calculate 
luminosity density error values. No extinction correction has 
been applied. 

In the redshift range z > 4 the data points show 
fairly large discrepanc i es. As has been pointed out in 
ISawicki and Thompsonl (|2006bl ). the value obtained by the 
SDF for z ~ 4 is about 2.5 times larger than that by the 



KDF for the integration down to L = 0.1L^ =3 and the dis- 
agreement reaches a factor of 3.7 when extrapolated to zero 
luminosity. The situation is similar at 2 ~ 5 where SDF 
estimates of UV luminosity density are larger than our re- 
sul ts by factors of 2-3. The large luminosity density values 
by lOuchi et~aH j|2004al ) come from relatively large number 
densities in the fainter part o f their UVLF (see Fig. 9 of 
ISawicki and Thompson|[2006al ). By unconstrained Schechter 
function fitting they obtained the best-fitted faint end slope 
of a — —2.2 fo r z ~ 4 LBGs, which is signific antly steeper 
than the one bv lSawicki and Thompsonl (|2006al ) (-1.26) and 
the values estima t ed for 2^3 LBGs (« —1.6 to —1.4; 
ISteidel et al. I ll999l : ISawicki and Thompsonl [2006al ). In addi- 
tion to UV luminosity density estimates based on this steep 
q = —2.2 slope, they also examined cases with a fixed UVLF 
slope value of —1.6 to calculat e UV luminos i ty den sities at 
2 = 4-5, and all results from lOuchi et al.l l|2004ah shown 
in Fig. [TS] are based on calculations assuming the faint end 
slope fixe d to —1.6. 

The lOuchi et all l|2004ah z ~ 5 number density at 
Muv ~ —20.8 mag (the faintest data point for them) is com- 
parable to ours, but since their brighter part of the UVLF 
has smaller number densities than does our LF, their best- 
fitted M* is fainter and 0* is larger than ours. Their assumed 
faint-end slope ( — 1.6) is also steeper than our best-fit value 
( — 1.48). As a result, the contribution from fainter objects 
(extrapolated beyond the observed limiting magnitude) to 
the UV luminosity becomes large r than that in our d ata. 
The large differences between the lOuchi et al.l l|2004al) to- 
tal UV luminosity densities and those integrated down to 
0.1L* =3 , at both 2 ~ 4 and 2 ~ 4.7, also indicate that in 
their results the UV luminosity density comes mainly from 
unexplored faint populations. Because the bulk of UV lu- 
minosity emitted into inter-galactic space is from galaxies 
wi th UV luminosity L uy 5- L* , luminosity ranges probed 
bv lOuchi etail (|2004aT) and |I03| are too shallow to definitely 
measure the total UV luminosity. This is also the case for 
our result presented in this paper, since our data are deeper 
but still only reach to 0.3L* =3 . For 2 ~ 6, the studies pre- 
sented here both (at least partly) use the common Hubble 
Ultra Deep Field (UDF) data but, despite this, their results 
differ from each other. To discuss in detail what causes this 
discrepancy at 2 ~ 6 is beyond the scope of this paper, 
but considering the facts that 2 ~ 6 LBGs are primarily 
selected by one i — 2 coloui0, and that in both studies UV 
continuum luminosity estimated from z'-band magnitudes 
may suffer attenuation by intergalactic neutral hydrogen, 
this diversity at z ~ 6 may be caused by the uncertainty 
in sample selection and /or cor rection for UV luminosity (cf. 
ISawicki and Thompsonll2006bl ). 

Despite the uncertainty in the higher redshift ranges, it 
would be meaningful to examine the overall change of UV 
luminosity density with redshift using the currently avail- 
able data in order to see the cosmic evolution of sta r forma - 
tion activity. The trend suggested bv lGabasch et ail (|2004bl ) 
based on their photometric redshift sample, is that at 2 > 3 
the UV luminosity density is slowly declining with increas- 
ing redshift. This trend is consistent with the combination of 



3 Some p art of UDF has very deep near-infrared imaging with 
NICMOS jBouwens et al.ll2006l) . 
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Figure 15. UV luminosity densities at various redshifts. Coloured symbols: UV luminosity densities integrated down to zero luminosity. 
Black-and-white symbols: UV luminosity densities integrated down to Lyy = 0.1L*_ 3 (Muy = —18.5 mag). Filled circles at z = 4.8: 
this work, using the UVLF derived from the combination of our two survey fields. Open triangles at < z sC 1: based on UVLFs by 
Takeuchi et al. (2005). Filled squares at 2.2 < z < 4: Sawicki and Thompson (2006b). Crosses at z = 4.0 and z = 4.7: Ouchi et al. 
(2004). Filled triangles at z = 6: Bouwens et al. (2006). Pentagon at z = 6: Bunker et al. (2006). Open circles: Gabasch et al. (2004). For 
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ranges (Ttjv ^ 0.1L*_ 3 and L > 0). The right vertical axis indicates the comoving star formation density converted from UV luminosity 
density. The conversion factor in the logarithmic scale is 27.903 (Madau et al. 1998). 



ISawicki and Thompson! (l2006a|) and our resu l ts, and the UV 
luminosity density values of iGabasch et al.l l|2004bl ). deter- 
mined from the UVLF for galaxies with photometric redshift 
4 < z < 5, also agree reasonably well with ours results. If 
we restrict ourselves to see only data points with integration 
down to 0.114=3 (black and white points in Fig. [15}, which 
are less affected by the uncertainty of the UVLF faint-end 
slope, then all the data in Fig. [15] show a consistent picture 
that the UV luminosity density (i.e., comoving star forma- 
tion density, if we neglect the effect of dust attenuation) 
is low at z ~ 6 (when the cosmic age is pb 0.9 Gyr) and 
gradually increases until z ~ 2-3 (when the cosmic age is 
2-3 Gyr). From that epoch on the UV luminosity density 
declines rather rapidly until now. 

This overall cosmic evolution of UV luminosity density 
might indicate that some kind of transition in star-forming 
galaxies has happened at z ~ 2 — 3. At z < 2, it has been 
pointed out that in the later epochs there is an absence of 
massive galaxies with large SFR while less massive galaxies 
continue to form stars. This trend is called "down-sizing" 
and recent development of deep and wide imaging surveys 
and intensive surveys using multi - objec t spectrographs (e.j 



spec 

Kodamagt all 12004 Ifreu et al.l 120051 : I Juneau et aD [200 



Perez-Gonzalez et al ] l2005l )_ made this trend clearer than 



when it was first suggested ijcowie et al J 199*3 ) . On the other 
hand, in this paper we have argued that there is differential 
galaxy evolution at z > 3, i.e., that the galaxy population 
with high SFRs maintains its number density from early 
times, z ~ 5, until z ~ 3, while the number density of galax- 
ies with smaller SFRs gradually increases. Moreover, we also 
discussed that this differential evolution could be an indica- 



tion of biased galaxy evolution in the early universe. The 
break in the UV luminosity density around z «2-3 may in- 
dicate the transition from biased galaxy evolution to down- 
sizing. In other words, the formation of massive galaxies, i.e., 
the active assembly of stellar mass within the massive dark 
matter haloes, might have started at z > 6 and continued 
at roughly constant rate for more than 2 Gyr, with the bulk 
of stellar mass assembly in such massive haloes terminating 
around z ~ 2. Although the observational evidence of down 
sizing has accumulated rapidly, the biased galaxy evolution 
at higher redshift is still just a hypothesis. Moreover, cur- 
rently we only deal with rest-frame UV properties which are 
dominated by the effects of short time-scale starbursts. In 
order to investigate what is happening at z > 3 and to relate 
that activity to down-sizing in the lower redshift range, we 
need information on masses (either stellar or of dark matter 
haloes) . 



6.3 Contribution to the Ionizing Background 
Radiation 

UV luminosity density can be used to estimate the contri- 
bution of stellar sources to the cosmic ionizing background 
radiation under an assumption about the ratio of ionizing 
radiation fl ux e scaping from galaxies to non-ionizing UV ra- 
diation. In Il03l we argued that if the escape fraction is as 
large as 0.5, then ionizing radiation from star-forming galax- 
ies and from QSOs at 2 ~ 5 may be sufficient to keep the 
universe at that redshift ionized. Several studies of galaxies 
at z ~ 6 have discussed whether the UV luminosity den- 
sity of z ~ 6 galaxies is sufficient to keep the universe 
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ionized (e.g., Bouwens et al. 20031; Bouwens et ah! 120061 : 



lYan and Windhorst! |2004| ; iBunker et all l2006h . However, 
since the escape fraction of ionizing photons from galaxies 
is still almost unknown observationally (e.g , Steide l et al 
200 ll ; iGiallongo et"aill2002l ; llnoue et al . 20051: IShaplev et al 



2006), such arguments cannot be definitive. Additionally, 



the uncertainty in the UVLF, especially in its faint-end 
slope, makes th e estimates of ioniz ing radiation difficult. On 
the other hand, llnoue et all l|2006r i used the evolution of the 
UV luminosity density presented in Fig. If 51 emission from 
QSOs, 1GM opacity model and ionizing background radia- 
tion in the redshift range to 6 to predict the escape frac- 
tion of star-forming galaxies. They found that evolution of 
the escape fraction — i.e., higher (> 0.2) escape fraction at 
higher redshift (z ^ 4) — may be required. Unfortunately, 
the direct observation of ionizing radiation at z > 4 is ex- 
tremely difficult due to the increasing opacity of the 1GM. 
We should first investigate the properties of ionizing radia- 
tion from galaxies at z < 3 to know the typical value of the 
escape fraction in the high-redshift universe and its relation 
to other physical properties of galaxies. 



7 SUMMARY AND CONCLUSIONS 

In this paper we have reported the construction of a large 
sample of Lyman break galaxies (LBGs) at z ~ 5 and their 
UV luminosity function (LF) , obtained from deep large-area 
Subaru observations in two blank fields, namely, the Hubble 
Deep Field-North and the J0053+I234 regions. Our main 
findings are as follows: 

• Over an effective survey area of 1290 arcmin 2 , we ob- 
tained 853 LBG candidates at z ~ 5 with z' < 26.5 mag. 
Spectroscopic redshift data for galaxies at lower redshifts 
(z < 4) from published redshift surveys in the target fields 
as well as the results of follow-up spectroscopy for a subset 
of our LBG candidates confirm that our colour selection cri- 
teria are defined so that star-forming galaxies at z ~ 5 are 
efficiently selected and contaminations from objects at lower 
redshift are properly eliminated. 

• In the bright end of the UVLF there is no significant 
change in the number density from z ~ 5 to 3, while we 
found a gradual increase from z ~ 5 to 3 in the LF's 
fainter part (Muv > —21 mag). Such differential evolu- 
tio n of the UVLF from z ~ 4 to 3 has been pointed out 
by ISawicki and Thompson! l|2006ah , and our present work 
shows that this trend was already underway at earlier cos- 
mic times. 

• We found a deficiency of galaxies with blue UV colours 
(I c - z' < 0) at bright magnitudes (z < 24.5 mag). This 
trend is significant even if we take photometric errors into 
account. 

• The UV luminosity density at z ~ 5, calculated by in- 
tegrating the UVLF, is found to be 38.81 \l% of that at 
z ~ 3, for the luminosity range L > 0.1L* =3 . By combining 
our results with those from the literature, we found that the 
cosmic UV luminosity density is slowly increasing with cos- 
mic age from z ~ 6, that it then marks its peak at z = 3-2, 
and then declines until the current epoch. 

We have discussed the origins of the differential evolu- 
tion of the UVLF along cosmic time and suggested that our 



observational findings are consistent with the biased galaxy 
evolution scenario. A galaxy population hosted by massive 
dark haloes starts active star formation preferentially at 
early times and maintains its activity until z ~ 3, while 
less massive galaxies increase their number density later. 
This idea is supported by other observational findings, in- 
cluding, (I) the deficiency of UV luminous objects with blue 
UV colours, (2) the fact that UV luminous objects at z ~ 5 
have smaller Lyman a equivalent widths compared to fainter 
objects, and (3) the dependence of clustering strengths on 
UV luminosity which suggests that UV luminous objects are 
more likely to be hosted by massive dark matter haloes. Fur- 
ther investigations of the properties of LBGs, especially at 
higher redshifts (z > 4) and with wider dynamic ranges in 
luminosity and wavelength, are required to test the reality 
of the biased galaxy evolution scenario. 
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Table 1. A summary of observations for HDF-N. Limiting magnitude is for 5 cr with 1.2"diameter aperture. 



Filter Single Exposure Time 


Number of shots 


Total integration time 


Limiting mag. 


Seeing 


sec 




sec 




// 


V 720-1,440 


14 


24,720 


28.17 


1.11 


600 


10 








Ic 180-300 


11 


6,180 


26.86 


1.13 


360-420 


10 








z' 120-215 


21 


11,935 


26.55 


1.13 


180-280 


41 








Table 2. A summary of observations for J0053+1234. Limiting magnitude 


is for 5 a with 


1.2"diamet 


aperture. 










Filter Single Exposure Time 


Number of shots 


Total integration time 


Limiting mag. 


Seeing 


sec 




sec 




// 


V 600 


23 


13,800 


27.79 


0.89 


I c 200 


10 


2,000 


26.39 


0.87 


z' 240 


28 


6,720 


26.24 


0.88 



Table 3. Cross identifications of objects with previously published spectroscopic redshifts at 4.5 < 
in the HDF-N region. 



< 5.5 



ID 


R.A.(J2000) a 


Dec.(J2000) a 


z' mag. 


V - 


Ic 


Ic- 


z< 


Redshift 


Ref. 6 


F36219-1516 


12:36:21.88 


+62:15:17.0 


25.39 





67 





09 


4.890 


1 


F36279-1750 


12:36:27.74 


+62:17:47.8 


26.05 


>2 


32 


-0 


11 


4.938 


1 


F36376-1453 


12:36:37.62 


+62:14:53.8 


21.97 


3 


25 





78 


4.886 


1 


HDF 4-439.0 


12:36:43.84 


+62:12:41.7 


25.29 


1 


48 


-0 


19 


4.54 


2 


HDF 4-625.0 


12:36:44.65 


+62:11:50.7 


25.24 


>2 


87 





03 


4.580 


3 


GOODS J123647.96+620941.7 


12:36:47.93 


+62:09:41.7 


24.32 


>3 


31 


-0 


10 


5.186 


4 


ES 1 


12:36:49.23 


+62:15:38.8 


25.48 


>2 


96 


-0 


16 


5.190 


5 


HDF 3-0951.0 


12:37:00.23 


+62:12:19.8 


24.69 


1 


82 





22 


5.34 


6 


095819(A04-4) 


12:37:05.68 


+62:07:43.3 


24.50 


2 


26 





15 


4.650 


7 


GOODS J123721. 03+621502.1 


12:37:21.00 


+62:15:02.1 


23.32 


2 


36 





56 


4.761 


8 


144200(A04-2) 


12:37:57.49 


+62:17:19.0 


24.08 


2 


40 





07 


4.695 


7 


129178(A04-7) 


12:38:04.36 


+62:14:19.7 


24.29 


>3 


28 





21 


5.183 


7 


104268(A04-1) 


12:38:11.32 


+62:09:19.4 


24.02 


2 


22 





11 


4.517 


7 


148198(A04-8) 


12:38:16.63 


+62:18:05.3 


24.50 


>3 


19 





12 


4.615 


7 


149472 (A04-6) 


12:38:25.52 


+62:18:19.7 


24.87 


>3 


02 





11 


4.857 


7 


139294(A04-5) 


12:38:28.96 


+62:16:18.8 


24.39 


1 


96 





20 


4.667 


7 



a Positions are based on our astrometry with Suprime-Cam images, using USNO-B1 catalog (see section 
3.2) 

b 1: iDawson et al.l ll200lh. 2: Istern and Spinradl il999l) . 3: iFernandez-Soto et all Jl999T). 4: iBarger et alj 
l|2002h . 5: lDawson et a.1.1 P002I V 6: ISpinrad et aiTl|l99Sft . 7: lAndo et al] j2004h . 8: ICowie et al] j2004h . 



Table 4. Cross identifications of objects with previously published spectroscopic redshifts at 4.5 < z < 5.5 
in the J0053+1234 region. 



ID 


R.A.(J2000) a 


Dec.(J2000) a 


z' mag. 


V -Ic 


Ic-Z' 


Redshift 


Ref. 6 


106426(A07-1) 


00:52:21.34 


+12:32:35.3 


24.07 


2.59 


0.23 


4.797 


1 


093014(A07-4) 


00:52:37.37 


+12:29:58.7 


25.26 


>1.85 


0.16 


4.491 


1 


091813(A07-3) 


00:52:39.88 


+12:29:44.1 


25.03 


1.98 


0.06 


4.391 


1 


104115(A07-2) 


00:52:43.27 


+12:32:08.2 


24.22 


1.75 


0.16 


4.267 


1 


CDFa G01 


00:53:33.21 


+12:32:07.3 


23.54 


2.88 


0.03 


4.815 


2 


CDFa GD07 


00:53:35.57 


+12:31:44.1 


23.72 


1.67 


0.17 


4.605 


2 


CDFb G05 


00:53:51.28 


+12:24:21.3 


24.31 


2.70 


0.29 


4.486 


2 



a Positions are based on our astrometry with Suprime-Cam images, using USNO-B1 catalog (see section 
3.2). 

6 l: |Ando et all 1120071 1. 2: ISteidel etHl dl999h . 
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Table 5. Number of LBG candidates and UV luminosity function of LBGs at z ~ 5 for HDF-N. 



z' 


mag 


abs. mag 


Num 


V efI 


N int 


0(m) 












Mpc~ 3 




N/Mpc _3 /mag 




23.0 


- 23.5 


-23.082 


2 


3.07 x 10 3 


0.79 


1.55 ±0.85 x 10" 


-6 


23.5 


- 24.0 


-22.582 


5 


2.86 x 10 3 


2.79 


3.04 ±0.91 x 10" 


-6 


24.0 


- 24.5 


-22.082 


35 


2.68 x 10 3 


6.06 


4.25 ±0.66 x 10" 


-5 


24.5 


- 25.0 


-21.582 


72 


2.48 x 10 3 


12.58 


9.42 ± 1.05 x 10" 


-5 


25.0 


- 25.5 


-21.082 


126 


2.24 x 10 3 


24.94 


1.78 ±0.15 x 10" 


-4 


25.5 


- 26.0 


-20.582 


186 


1.90 x 10 3 


24.05 


3.36 ±0.25 x 10" 


-4 


26.0 


- 26.5 


-20.082 


191 


1.51 x 10 3 


20.33 


4.43 ±0.33 x 10" 


-4 



Table 6. Number of LBG candidates and UV luminosity function of LBGs at z ~ 5 for J0053±1234. 



z' 


mag 


abs. mag 


Num 




N int 


0(m) 












Mpc~ 3 




N/Mpc _3 /mag 




23.0 


- 23.5 


-23.082 


3 


2.51 x 10 3 


1.40 


1.63 ±0.69 x 10" 


-6 


23.5 


- 24.0 


-22.582 


10 


2.38 x 10 3 


1.08 


9.61 ± 2.88 x 10" 


-5 


24.0 


- 24.5 


-22.082 


22 


2.28 x 10 3 


2.53 


2.19 ±0.44 x 10" 


-5 


24.5 


- 25.0 


-21.582 


79 


1.99 x 10 3 


7.82 


9.16 ± 1.02 x 10" 


-5 


25.0 


- 25.5 


-21.082 


122 


1.68 x 10 3 


14.85 


1.63 ±0.15 x 10" 


-4 



Table 7. UV luminosity function of LBGs at z ~ 5 obtained as a weighted average of results for two survey 
fields. The number densities at z' > 25.5 come only from the HDF-N region. 



z 1 mag abs. mag <l>{ m ) 

N /Mpc - 3 /mag 



23.0 


- 23.5 


-23.082 


1.60±0.75 x 10" 


-6 


23.5 


- 24.0 


-22.582 


7.02±2.11 x 10" 


-6 


24.0 


- 24.5 


-22.082 


3.00±0.53 x 10" 


-5 


24.5 


- 25.0 


-21.582 


9.26±1.03 x 10" 


-5 


25.0 


- 25.5 


-21.082 


1.69±0.15 x 10" 


-4 


25.5 


- 26.0 


-20.582 


3.36±0.25 x 10" 


-4 


26.0 


- 26.5 


-20.082 


4.43±0.33 x 10" 


-4 



Table 8. Schcchtcr function parameters of UVLF at z ~ 5, obtained with the data from the HDF-N region 
only and from both the HDF-N and J0053+1234 regions. 



Region 


M* 


4> 


* 




a 


HDF-N 


-20.82t°-g 


7 7+1.9 
'■'-2.1 


X 


io- 4 


-1 05+ - 35 
1 - UO _0.32 


HDF-N and J0053±1234 


-2l.28±°; 33 


4 1+2-9 


x 


10- 4 


4 o+0.38 
± -^°-0.32 



Table 9. Evolution of the number density for galaxies brighter and fainter than M\jy = —21.0. 



Epochs 


bright end (M uv < -21.0) 


faint end (M uv > -21.0) 




^(higher z) /^(lowerz) 


^(higher z) /^(lowerz) 


2 ~ 5 — > 3 


1.37 ± 0.44 


0.34 ± 0.02 


z ~ 5 — > 4 


1.58 ± 0.56 


0.78 ± 0.09 


z ~ 4 -> 3 


0.87 ± 0.19 


0.44 ± 0.05 
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Table 10. UV luminosity density calculated for the luminosity range Ltjv ?J 0.1L*_ 3 (Mtjv —18.5) and 
L uv > 0. 



redshift 


L vv > 0.1L* =3 


iuv > 


Reference 


5 n+ - 5 

°' U -0.5 

5 n+ - 5 

°- U -0.5 


25 81 + ' 04 
^o.»i_ 04 

95 84+ - 06 


25.86±°;£ 


HDFN 

HDFN and J0053+1234 



